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Chapter  1 
INTRODUCTION 


A  useful  tool  for  the  designer  of  an  optical  system  is  a 
technique  to  determine  the  effect  of  deviations  from  the  design 
positions  of  both  the  input  beam  and  the  system's  optical  surfaces. 
This  report  develops  a  technique  for  tracing  a  light  beam  through 
an  optical  system  for  application  on  a  digital  computer.  The 
optical  system  is  comprised  of  a  large  class  of  reflecting  and 
refracting  surfaces  that  include  cylindrical  optics.  The  ray 
tracing  technique  is  restricted  to  geometrical  optics  and  does  not 
require  rotational  symmetry  to  specify  a  quadric  surface. 

RAYTRACE  is  the  name  of  the  computer  program  that  traces  the 
light  beam  through  the  optical  system.  As  a  light  beam  is  traced 
through  an  optical  surface ,  the  number  of  rays  that  have  failed 
to  pass  the  surface  is  given.  The  rays  that  are  outside  a  given 
minimum  and  maximum  radius  at  the  output  distance  from  the 
surface  are  considered  to  be  ray  failures.  After  the  beam  ha3 
been  traced  through  the  optical  system,  the  beam  can  be  retraced 
through  the  system  with  new  input  data  by  entering  the  data 
through  the  console  typewriter.  The  output  options  include  a 
listing  of  the  output  beam  parameters  and  graphics  at  the  output 
of  any  particular  optical  surface.  The  graphics  consist  of  a 
cross  sectional  view  and/or  knife  edge  scan  of  the  output  beam  at 
the  output  distance  from  the  optical  surface. 

A  plane  surface  with  the  optical  property  of  refraction  or 
reflection  is  the  basic  component  of  an  optical  system.  The 
relationship  between  the  ray  incident  to  the  plane  surface  and 
the  ray  refracted  or  reflected  from  the  plane  surface  is 
developed  for  a  plane  surface  system.  This  relationship  is 
referred  to  as  t'ie  input-output  ray  equations.  Since  most  mirrors 
and  lenses  that  comprise  an  optical  system  are  quadric  surfaces 
rather  than  planar,  the  relationship  between  a  quadric  surface 
and  a  plane  surface  must  also  be  established.  Definitions  of  the 
parameters  and  symbols  used  throughout  the  text  are  listed  in  the 
^lousary,  Appendix  A. 
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MATHEMATICAL  BACKGROUND 

Tha  coordinate  systems  defined  in  this  report  are  right 
handed  orthogonal  coordinate  systems  whose  axes  are  Xj,  Y^,  and 

Zj.  The  subscript  j  designates  a  particular  coordinate  system 

where  j  *  1,  2,  ...  n.  An  alphabetic  subscript  in  conjunction 
with  a  numerical  subscript  designates  a  particular  point  in  a 
particular  coordinate  system.  For  example ,  point  A  in  the  1 
coordinate-  system  is  expressed  as  X^,  Y^A,  and  Z^.  Unit 

vectors,  which  are  parallel  to  the  X. ,  Y . ,  and  Z.  axes,  are 

A  A  A  j  j  j 

expressed  by  Xj,  Yj,  and  Zj.  Transformations  between  the  various 

coordinate  systems  are  determined  by  symbolic  representations  or 
Piograms  (ref.  (a)). 

The  direction  of  a  vector 

^  A  A  A 

V  *  AX  +  BY  +  CZ  (1) 

is  defined  by  the  three  angles  ct,  0,  and  y  measured  from  the 
positive  X,  Y,  and  Z  axes,  respectively,  to  the  vector  V.  The 
cosines  of  o,  0,  and  y  are  known  as  the  direction  cosines  of  the 
vector  V.  The  unit  vector  along  the  line  of  action  of  V  is 
given  by 

A  A  A  A 

V  *  cos  a  X  +  cos  0  Y  +  cos  y  Z.  (2) 

Since  the  magnitude  of  the  unit  vector  is  equal  to  one, 
obviously 

AAA 

cos*  a  +  cos*  0  +  cos  y  =  1.  (3) 

The  direction  of  a  vector  can  be  determined  by  the 
coordinates  of  two  points  (XA,  Yft,  Zft)  and  (x^,  Y0,  Zfi) .  The 

relationship  of  the  distance  between  the  two  points  and  the 
direction  cosines  can  be  expressed  by 

*B-  .  YB  -  YA  .  ZA  ... 

cos  a  cos  0  cos  y  '  1 

where  cos  a,  cos  0,  and  cos  y  are  not  equal  to  zero.  If  a 
denominator  is  zero  then  the  corresponding  numerator  is 
considered  to  be  zero. 

OPTICAL  BACKGROUND 

Since  an  optical  surface  can  be  either  refracting  or 
reflecting,  one  relationship  between  the  input  ray  and  the  out¬ 
put  ray  is  developed  to  handle  both  cases.  Straight  line 
propagation  is  assumed.  Taking  the  refraction  case  first,  as 
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shown  in  Figure  1,  the  refracting  plane  surface  is  designated 
by  the  Z  *•  0  plane,  and  the  incident  ray  has  direction 

A  AAA 

Vj  ■  cos  oij  X  +  cos  y  +  cos  Yj  Z  (5) 

where  the  subscript  I  denotes  the  incoming  or  incident  ray.  The 
refracted  ray  has  direction 

A  AAA 

VR  -  cos  oR  X  +  cos  0R  y  +  cos  YR  Z  (6) 

where  the  outgoing  or  refracted  ray  is  denoted  by  the  subscript  R. 

From  Figure  1  and  the  dot  product,  the  cosine  of  the  incident 
angle  is  given  by 

A  A 

COS  Tj  *  Z  •  Vj  *  -cos  Yj  (7) 

and 


sin  Tj  =  Vl  -  cos2  Tj.  «  fl  -  cos2  Yj*  (8) 

krom  Snell's  Law 


nT 

sm  t_  »  — —  sin  t_  (9) 

R  nR  I 

where  nj  is  the  index  of  refraction  of  the  medium  of  the  incident 
ray  and  nR  is  the  index  of  refraction  of  the  medium  of  the 
refracted  ray.  The  ratio  of  nx  to  nR  is  defined  by  the  symbol  N. 
From  equations  (8)  and  (9) 

cos  tr  *  Vl  -  sin2  tr  *  Vl  -  N2  (1  -  cos2  Yj)  •  (10) 

Since  the  direction  cosine  is  measured  from  the  positive  axis 

cos  tr  »  cos  (180°  -  tr)  *  -cos  tr  (11) 


where  cos  yr  is  the  direction  cosine  of  the  refracted  ray. 
Therefore, 

cos  yr  =  -  Vl  -  N2  (1-cos2  Yj)* 


(12) 
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Since  the  component  of  the  ray  parallel  to  the  refracting 
surface  does  not  change  direction,  it  must  be  true  that 


From  equation  (3) 


cos  aR  cos  otj 
cos  'Br  *  cos'Bj* 


cos2  aR  +  cos2  BR  +  cos2  yR  *  1. 


(13) 


(14) 


Using,-  equations  (12)  and  (14) 


cos2  aR  +  cos2  BR  -  1  -cos2  yR  *  N2  (1-coo2  yT). 


Utilizing  equation  (13) 


2 

cos  Bt 


(cos2  Oj \ 

zzrz) 


+  cos2  $R  *  N2  (1-cos2  y^) 


2 

cos  Bt 


( 


cos2  Oj  +  cos2  Bj 
- TZ - 

COS  Bt 


^  =  N2  (1— cos2  Yj) 


cos2  B 


ii 


f  _  :°32  m  * 

\  cos^  Bt  / 


(1-cos  Yx) . 


Therefore , 


or 


By  the  same  method 


cos2  Br  *  N2  cos2  Bj 


cos  BR  *  N  cos  Bj. 


cos  aR  »  N  cos  Oj. 


Therefore,  the  direction  of  the  refracting  ray  is 


(15) 


VR  “  C03  <XR  X  +  COS  Bp  Y  +  COS  Yp  Z 


(16) 
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where 


cos  aR  *  N  cos  cij  ( 

cos  8r  *  N  cos  Bj  ( 


cos  yR  *  -  f 1  -  N2  (1-cos2  Yj) .  ( 

For  a  reflective  surface,  the  angle  of  the  incident  ray  is 
equal  to  the  angle  of  the  reflected  ray  as  shown  in  Figure  2.  ' 

obtain  the  reflected  ray  direction,  the  sign  of  the  normal 
component  of  the  incident  ray  in  equation  (5)  i3  reversed. 
Therefore,  the  relationship  of  an  incident  ray  with  direction 

A  AAA 

Vj  =*  cos  ttj  X  +  cos  Y  +  cos  Yj  Z 
and  a  refracted  or  reflected  ray  with  direction 

/V  AAA 

VR  *  cos  aR  X  +  cos  BR  Y  +  cos  yr  Z 
is  obtained  by  rewriting  equations  (17) 

cos  aR  *  N  cos  aR  C 


where 


cos  Br  *  N  cos  8r 


cos  yr  -  -RP 


V 1  -  N2  (1-cos2  yt) 


RP 


1 

-1 


refraction 

reflection 


N 


nj/nR  refraction 
1  reflection 


In  developing  the  equations  establishing  the  relationship 
between  the  input  rays  and  the  output  rays,  the  complement  to 
the  angle  of  incidence  and  the  complement  to  the  angle  of 
refraction  or  reflection  is  used.  If  Bj  is  the  complement  of 

the  angle  of  incidence,  then  BR,  the  complement  of  the  angle  of 

refraction,  can  be  determined  from  Snell's  Law. 


sin  (tt/2-8-)  =  N  sin  (it/2-Bv) 
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or 


cos  ©R  »  N  cos  ©j. 


Therefore,  ©K  for  the  output  ray  can  be  determined  by 
«* 

equation  (19)  as  follows 

©R  *  RP  cos"*3  (W  cos  ©j) 


where 


|N  COS  ©jl  <  j 

for  the  refracted  ray  to  exist • 


(19) 

rewriting 

(20) 
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Chapter  2 

PLANE  SURFACE  SYSTEM 


REFERENCE  COORDINATE  SYSTEMS 

The  plane  surface,  whether  it  be  reflecting  or  refracting, 
is  used  as  tne  basic  element  of  an  optical  system  in  this  report. 
This  plane  surface  along  with  the  associated  references  comprise 
the  plane  surface  system.  This  system  consists  of  a  reference 
input  ray,  a  reference  plane,  and  a  reference  output  ray.  The 
output  ray  is  either  a  refracted  ray  or  a  reflected  ray. 

Coordinate  system  1  is  the  coordinate  system  of  the  plane 
surface  reference.  As  shown  in  Figure  3,  the  plane  z^  =  0  is 

the  refracting  or  reflecting  surface  of  the  reference  plane  with 
the  origin  as  the  reference  point  of  incidence.  The  plane  of 
incidence  for  the  reference  input  ray  is  the  X^^  »  0  plane. 

Originally,  the  coordinate  system  attached  to  the  plane 
surface  is  coincident  with  coordinate  system  1.  However,  the 
plane  surface  can  move  both  in  translation  and  rotation  with 
respect  to  the  plane  surface  reference  system.  The  rotation 
angles,  e  and  t)»  are  measured  with  respect  to  the  coordinate 
system  attached  to  the  plane  surface.  The  coordinate  system 
attached  to  the  plane  surface  is  designated  coordinate  system  3. 
The  rotation  angles,  e  and  n*  are  defined  by  the  Piogram 
(ref.  (b)  and  »c))  of  Figure  4.  The  coordinate  (X1D,  Y1D,  Z1D) 

describes  the  translation  of  the  origin  of  the  plane  surface  from 
its  reference  position,  the  origin  of  coordinate  system  1.  The 
Z3  *  0  ylane  is  the  refracting  or  reflecting  surface  after  the 
translation  and  rotation  of  the  plane  surface. 

The  coordinate  system  of  the  reference  input  ray  as  shown  in 
Figure  3  is  coordinate  system  2.  The  origin  of  this  coordinate 
system  is  located  at  a  distance  from  the  origin  of  system  1 

along  the  reference  input  ray  with  direction  Zy  The  =  0 

plane  is  coincident  with  the-  X2  3  0  plane. 

Coordinate  system  4  is  the  coordinate  system  of  the  reference 
output  ray.  The  origin  of  coordinate  system  1  and  }.he  origin  of 
coordinate  system  4  are  separated  by  the  distance  RR  along  the 

output  reference  ray  with  direction  Z^.  The  =  0  plane  is 

coincident  with  the  =  0  plane. 
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Input  Reference  Plane 


Figure  3  Reference  Planes  for  Plane  Surface  /  /  Y,  \  v 
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Input  Ray  Deviation 

The  input  ray  deviation  from  the  input  r&hrrence  ray  with 
direction  Z„  can  be  uniquely  described  by  two  ).inear  positions 

d* 

and  two  angular  positions.  The  linear  position  is  simply  the 
distance  from  the  origin  in  the  Z2  *  0  plane.  ?he  coordinates 

of  the  input  ray  deviation  from  the  input  reference  ray  are 
(X2D#  Y2d)  .  The  angular  deviation  from  the  inj>5  .t  reference 

direction  is  the  deviation  from  the  Z2  direction.  These  two 

angular  deviations  from  the  input  reference  direction  are  02 

and  4*2*  As  shown  in  Figure  5  02  is  the  angle  that  the  projection 

of  the  input  ray  on  *  0  plane  makes  with  the  Z2  direction. 

4»2  is  the  angle  that  the  projection  of  the  input  ray  on  the  X2 

*  0  plane  makes  with  the  Z2  direction. 

The  direction  cosines  of  the  input  ray  can  be  determined 
from  the  input  ray's  deviation  angles  and  equation  (3).  From 
Figure  5 

X  *  V  cos  ct2  (21) 

Y  =  V  cos  B2 

Z  *  V  cos  y2 

and 

tan  02  =  X/Z  (22) 

tan  2  =  Y/z* 

Therefore#  by  substitution  of  equation  (21)  in  equation  (22) 


tan  02  -  cos  a2/cos  Y2» 
tan  \j<2  =  cos  B2/cos 

Equation  (3)  can  be  written  as 

2  2  2 
cos  a2  +  cos  B2  *  sin  y2 


2  2  .2 

cos  a2  cos  B2  sm  If 2 

- 5 -  +  - 5 -  =  - 2 - 

cos  y2  cos  y2  cos  >2 


(23) 


(24) 
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Using  equation  (23)  in  equation  (24) 

tan2  02  +  tan2  ij>2  =  tan2  y2«  (25) 

Hence,  equations  (23)  and  (25)  can  be  used  to  determine  the 
direction  cosines  of  the  input  ray.  The  same  considerations 
apply  to  the  output  ray. 

Input-Output  Ray  Equations 


The  linear  deviations  of  the  input  ray  from  the  input  reference 
ray  are  X2D  and  Y2D.  The  input  ray  direction  is 

A  A  A  A 

Vj  =•  cos  a2  X2  +  cos  B2  Y2  +  cos  y2  Z2  (26) 

with  angular  deviations  of  j?2  and  ^2.  The  plane  surface  translations 
are  X1D,  Y1D,  and  Z1D  with  angular  rotations  of  e  and  n. 


The  Piogram  of  Figure  6  relates  the  input  reference  system 
(coordinate  system  2)  to  the  coordinate  system  attached  to  the 
plane  surface  (coordinate  system  3) .  Since  unit  vectors  in 
rectangular  coordinate  systems  are  independent  of  position,  the 
Piogram  of  Figure  '5  is  used  to  determine  the  relationships  emong 
the  unit  vectors  by  setting  Rj  =  X1D  =  Y1Q  *  Z1D  »  0.  Then  the 

input  ray  direction  with  respect  to  coordinate  system  3  is 
determined  from  these  relationships. 

The  transformation  of  vectors  from  system  2  to  system  3  is 
given  by 

'cos  n  cos  (Oj+e)  sin  n  sin  (6j+e)  sin  n 
Tj  =  0  -sin  (Gj+e)  cos  (Sj+e)  (27) 

„sin  n  -cos  (Gj+e)  cos  n  -sin  (G^+e)  cos  r)- 


Using  equation  (27) ,  the  direction  of  the  input  ray  with  respect 
to  system  3  is 

A  A  A  A 

Vj  *  cos  «3  X3  +  cos  S3  Y3  +  cos  y3  Z3  (28) 

where 
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The  coordinate  of  the  point  where  the.  input  ray  intersects 
the  input  refersr.ee  plane  in  terms  of  syr.tem  .>  coordinates  is 
denoted  by  the  subscript  Q.  The  coordinate  relationship  can 
also  be  determined  from  the  Piogram  in  Figure  6  with  the  equation 
given  as  follows: 


X3Q 

X2D 

cos  ri  sin  e  sin  n  -cos  e  sinn 

X1D 

Y3Q 

»  T, 

- 

0  cos  sin  e 

yid 

.Z3Q. 

-Rx 

«■ 

sin  n  -sin  c  cos  n  cos  e  cos  n 

.Z1D. 

Using  equation  (4)  the  equations  for  the  point  of  incidence 
on  the  plane  surface  in  system  3  can  be  determined.  The  point  of 
incidence  is  denoted  by  the  subscript  P  and  the  plane  of  incidence 
is  given  by  the  plane  Z,  =  0.  Therefore,  the  point  of  incidence 
on  the  plane  surface  is'' 


X3P  -  * 


Y3P  -  1 


3Q  “ 

3Q  cos 

3 

Y3 

Ola) 

cos 

2  _ 

*3 

(31b) 

3Q  “ 

3Q  cos 

Y3 

Z3P  *  0 


In  order  to  determine  the  direction  of  the  output  ray 
equation  (18)  is  used.  In  addition,  equation  (20)  is  used  to 
determine  the  complement  of  the  angle  of  refraction  or  reflection. 
The  subscript  R  denotes  the  output  ray  and  the  prime  (’)  refers 
to  the  output  ray  direction.  The  direction  of  the  output  ray 
with  respect  to  system  3  is 

A  ^  A  A 

VR  *»  cos  a*  3  x3  +  cos  8*2  Y3  +  cos  y*  3  Z3  (32) 


where 


cos  a*3  =  II  cos  a3 
cos  8'3  =  II  cos  33 

cos  y  '  2  =  -RP  Vl  -  i;2  (1-cos2  y3) 


and 


6 


R 


RP  cos 


CJ  cos  Oj) 


(30) 
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In  Figure  7  the  Piogram  shows  the  relationship  between  the 
coordinate  system  attached  to  the  plane  surface  -{system  3)  and 
the  output  reference  plane  (system  4).  The  transformation  of 
vectors  from  system  3  to  system  4  is 


cos  n 


sin  n 


Tr  *  sin  n  cos  (SR+e)  -sin  («R+e)  -cos  n  cos  («R+e) 

Lain  n  sin  (6R+e)  cos  (6r+e)  -cos  n  sin  (©R+e)_ 

The  direction  of  the  output  ray  with  respect  to  system  4  is 


cos  a4  X4  +  cos  8^  VA  +  cos  yA  Z 


4  *4 


4  “4 


where 


cos  a. 


cos  8, 


cos  Y; 


cos  a' 


cos  8*3 
cos  y'3 


The  coordinates  of  the  point  where  the  output  ray  intersects 
the  plane  surface  in  terms  of  system  4  coordinates  are  given  by 


Y4Q  'TR  Y3P  '  0  sin9R  COs9R  Y1D  -  P  < 

LZ4oJ  LZ3fJ  1°  -°0S  9R  Sin  9R  J  LZ1dJ  LRrJ 

The  point  of  incidence  of  the  output  ray  on  the  output 
reference  plane  can  be  determined  by  using  equation  (4)  and  using 
the  fact  that  the  plane  of  incidence  is  given  by  the  plane  Z^  *  0 

Therefore,  the  point  of  incidence  of  the  output  ray  on  the  output 
reference  plane  is 

cos  a. 

X4D  ■  X4Q  '  Z4Q  (36a 


Y  -  V  _  r- 

4D  *4Q  *40 


Z4D  =  0 


cos  3^ 
cos  Yi 


The  output  ray  can  then  be  uniquely  described  by  two  linear 
positions  and  two  angular  positions  with  respect  to  the  output 
reference  plane  using  equations  (23)  and  (36) .  The  linear  and 
angular  positions  are 
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cos  a. 

a4D  a4Q  4Q  cos  y4 

(37a) 

cos  6. 

v  —  v  _  n  ** 

4D  *4Q  *4Q  cos  Y4 

(37b) 

04  =  tan-1  (cos  o4/cos  y4) 

(37c) 

'J>4  =  tan-1  (cos  84/cos  y4>  • 

(37d) 

These  output  ray  parameters,  X4D,  Y4Q,  0  4,  and  \J>4  are  used  as  input 
parameters  for  the  optical  surface  that  follows. 

INTERFACE  BETWEEN  TWO  PLANE  SURFACES 


When  two  plane  surfaces  are  in  series,  the  output  coordinate 
system  of  the  first  surface  may  not  be  coincident  with  the  input 
coordinate  system  of  the  second  surface  as  shown  in  Figure  8.  If 
the  two  coordinate  systems  are  not  coincident  an  interface  between 
the  two  systems  is  required.  This  interface  from  the  output 
coordinate  system  of  the  first  surface  to  the  input  coordinate 
system  of  the  second  surface  can  be  accomplished  by  a  rotation 
about  the  reference  ray  or  the  Z4  axis  as  shown  in  Figure  8.  A 

positive  rotation  is  determined  by  the  righthand  rule  about  the 
reference  ray  direction.  The  magnitude  of  the  rotation  is  the 
angle,  X,  between  the  positive  X  axis  of  the  output  plane  of  the 
first  surface  and  the  positive  X  axis  of  the  input  of  the  second 
surface.  The  Piogram  in  Figure  8  shows  the  relationship  between 
the  output  ray  parameters  of  the  first  surface  with  the  input 
reference  plane  of  the  second  surface. 


\ 


N 
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Chapter  3 
QUADRIC  SURFACES 


The  elements  that  comprise  an  optical  system  may  have  quadric 
surfaces  rather  than  plane  surfaces.  Therefore,  if  a  plane  surface 
is  considered  to  be  the  basic  subsystem  of  an  optical  system,  then 
the  relationship  between  a  quadric  surface  and  the  plane  surface 
system  must  be  developed.  The  basic  equation  of  a  quadric  surface 
can  be  written  as 

AX2  +  BY2  +  CZ2  +  DZ  =  E  (38) 

where  the  Z-axis  is  the  axis  of  symmetry.  The  center  of  symmetry 
may  translate  in  any  direction  with  respect  to  a  reference. 

However,  here  misalignment  of  the  axes  is  not  allowed.  The 
translation  of  the  center  of  symmetry  in  the  XY  plane  is  known  as 
decentration.  Focusing  is  the  translation  of  the  center  of  symmetry 
along  the  Z-axis.  Therefore,  equation  (38)  becomes 

A(X-C  ) 2  +  E  (Y-C  ) 2  +  C(Z  -  (C  +F))2  ♦  D(Z  -  (C  +F) )  -  E  (35) 
x  y  z  z 

where  C  and  C ,  are  the  decentration  coordinates ,  C  is  the  Z 
x  y  z 

coordinate  for  the  center  of  symmetry,  and  F  is  the  focusing  parameter. 

The  types  of  quadr ic  surfaces  generally  considered  for  an  optical 
system  are  an  ellipsoid ,  a  hyperboloid  of  two  sheets,  and  an 
elliptic  paraboloid  as  shown  in  Figure  9.  lhe  equations  for  these 
surfaces  are 


,2  2  2 

*  +  =  l  (ellipsoid) 

ar  b*  c 


(40  a) 


2  2  2 

-£j-^2  +  £j=1  (hyperboloid  of  two  sheets)  (40b) 
a  b  c 


2  2 

±  -  =  o  (elliptic  parabolcid) 
ar  d  c 


(40c) 


with  the  Z-axis  as  the  axis  of  symmetry.  Spheres  and  cylinders 
are  simply  special  cases  of  these  quadric  surfaces.  The  constants 
a,  b,  and  c  are  the  semi-axes  lengths  for  an  ellipsoid.  If  a  =  b  * 
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a)  Ellipsoid 


b)  Hyperboloid  of  Two  Sheets 


Figure  9  Quadric  Surfaces 
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a  sphere  is  the  result.  If  a  =  •»,  the  result  is  an  elliptic 
cylinder  with  elements  parallel  to  the  X-axis.  For  a  hyperboloid 
.of  two  sheets,  the  constant  c  is  the  vertex  of  the  surface  on  the 
Z-axis.  If  f  is  the  focus  of  the  hyperbola  in  the  Y  =  0  plane 

and  fy  is  the  focus  in  the  X  =  0  plane  then 

2  .2  2 
a  -  £x  -  c 


K2  =  f2  „2 

b  -  fy  -  C  . 

In  the  case  of  an  elliptic  paraboloid  the  constant  c  is  the 
location  of  an  ellipse  along  the  Z  axis  with  the  semi-axes  lengths 
of  a  and  b.  Table  i  lists  the  constants  that  are  substituted  into 
equation  (39)  to  obtain  the  desired  quadric  surface. 

A  relationship  between  a  quadric  surface  and  a  plane  surface 
system  can  be  developed  by  solving  for  the  linear  translation  and 
angular  rotation  of  the  plane  surface  from  the  reference  position 
to  a  point  tangent  to  the  input  ray  incident  on  the  quadric  surface. 
The  translation  and  rotation  parameters  are  equivalent  to  the 
coordinate  (X^D,  Y^D,  21D^  an<*  ^e  an9^es'  e  an<*  as  outlined  in 

Chapter  2.  Figure  10  shows  a  quadric  surface  with  the  reference 
coordinate  system  (system  1)  at  the  vertex  of  the  surface  and  an 
input  reference  system  (system  2) .  These  two  reference  systems 
are  separated  by  the  distance  R.  The  center  of  symmetry  for  the 
quadric  surface  is  the  Z  «  Cz  plane. 

Referring  to  Figure  10,  the  direction  of  the  input  ray  is 

A  A  A  A 

Vj  =  cos  <*2X2  +  cos  +  cos  y2Z2  (41) 

with  a  deviation  coordinate  of  (X2D,  Y2d' "  From  equation  (29)  and 
setting  6  =  0  =  0  and  0].  =  tt/2,  the  direction  cosines  and  the 
deviation  coordinate  with  respect  to  system  1  are 

cos  =  cos  c*2  (42) 


cos  3^  =  -cos  S2 

COS  -COS  Y2 


X1Q  =  X2D 
Y1Q  =  “Y2D 
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The  point  of  incidence  is  determined  from  equation  (4)  and  is 

X  *  +  (Z-R)  cos  cij/cos  (4 

Y  *  Y^q  +  (Z-R)  cos  Bj/cos 

where  (X,  Y,  z)  is  the  incident  coordinate  of  the  input  ray  on 
the  quadric  surface.  By  substituting  equation  (44)  into 
equation  (39) ,  the  result  is 

A  (G  +Z  tan  0-R  tan  0)2  +  B  (G  +Z  tan  lp-R  tan  \p)  ^  (4 

x  y 


+  C  (Z-G  ) *  +  D  (Z-G  )  *  E 

Z  Z 


where 


gy  =  Yio  "  cy 

Gz  =  cz  +  p 


tan  0  =  cos  a^/ccs 


tan  4<  -  cos  Bj^/cos  y^. 

By  combining  terms,  equation  (4S)  is  rewritten  as 

A?,2  +  %Z  +  Q  =o 


where 


A  »  A  tan2  0  +  B2  tan2  \j»  +  C 


6  =  2A  (G  -R  tan  0)  tan  0  +  2B  (G  -R  tan  \\>)  tan  ij;  -  2CG_  +  D 
a  y  z 

C  *  A  (G  -R  tan  0) 2  +  B  (G  -R  tan  ^)2  +  CG2  -  DG  -  E 
x  y  z  z 


with  two  conditions  that 


**  -  4A  >0 

1*4  |  +  18.  *  0 


(47a) 


(47b) 


If  either  of  these  conditions  are  not  met,  then  the  input  ray 
does  not  intersect  the  quadric  surface. 
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From  equation  (46) ,  the  z-coordinate  of  incidence  is 


Z 


+S  Vs2  -  AM 
2  A 

-?/f> 


At  t  0  and  B2  -  4Jip  >  0 
«&=  0  and  0 


(48) 


where  S  equals  ±1.  S  is  the  parameter  that  determines  whether  the 
quadric  surface  is  concave  or  convex.  In  order  to  solve  for  the 
correct  sign  of  Sr  let  0  «  i|»  =  F  =  Cx  *  Cy  =  X1Q  “  Z1Q  *  0. 

Substituting  into  equation  (48)  the  result  is 


2CG. 


-  D  +  S  Vd2  + 


4  EC 


2C 


G 

Z  D 


D*  +  4  EC  >  0  and  C  +  0 
C  =  0  and  D  ^  0 


(49) 


Using  Table  1  and  equation  (49)  the  result  for  an  ellipsoid  is 

Z  *  G  +  Sc. 
z 

Setting  the  center  of  symmetry  to  zero,  Z  must  be  positive  for  a 
convex  surface  and  negative  for  a  concave  surface.  Therefore, 

S  «  +1  for  a  convex  surface  and  S  =  -1  for  a  concave  surface. 

For  a  hyperboloid  of  two  sheets  the  equation  is 

Z  =  G  +  Sc 
z 

By  setting  the  center  of  symmetry  to  zero  the  surface  is 
concave  for  S  *  +1  and  convex  for  S  =  -1.  The  resulting  equation 
for  the  Z-coordinate  of  the  incident  ray  on  the  surface  of  an 
elliptic  paraboloid  is 


This  equation  provides  no  information  on  the  concavity  or 
convexity  of  the  surface.  This  information  can  come  directly 
from  equation  (40c).  If  c  <  0,  the  surface  is  concave  to  the 
input  ray  and  if  c  '  0  the  surface  is  convex  to  the  input  ray. 

A  summation  of  the  preceding  results  on  the  concavity  or 
convexity  of  a  quadric  surface  is  given  in  Table  2. 

Substituting  the  result  of  equation  (48)  into  equation  (44) 
yields  the  X  and  Y  coordinate  of  the  point  of  incidence  of  the 
input  ray  on  the  quadric  surface.  Now,  the  normal  to  the  surface 

is 

N  *  V  (A (X-C  ) 2  +  B (Y-C  ) 2  +  C (Z-G  ) 2  +  D (Z-G  )  -  E)  (50) 
x  y  z  z 
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and  the  aquation  for  the  unit  normal  is 

H  «  i  t(2A(X-C  ) }  X  +  {2B  (Y-C  )  }  Y  +  {2C(Z-G  )  +  D>  Z)  (51) 
«  a  y  z 

where 


N  *  f  4A2  (X-C  )2  +  4B2  (Y-C„)2  +  [2C  (Z-GJ  +  D]2.  (52) 

x  y  z 

Prom  the  Piogram  of  Figure  4,  the  normal  to  the  plane 
surface  is 

A  A  A 

=  sin  n  XL  -  sin  e  cos  n  Y1  +  cos  e  cos  n  Z^ .  (53) 

A  A 

Since  N  and  Z 3 are  both  normal  vectors  to  the  plane  surface,  the 
respective  components  of  the  unit  vectors  are  equivalent. 
Therefore, 


sin  n  “  2A  (X-Cx)/K 


or 

n  -  sin”1  (2A  (X-C  )/H) •  (54) 

Also, 

sin  e  cos  n  *  -2B  (Y-C  ,) 

x 

cos  e  cos  n  -  2C  (Z-G  +  D 
and  the  above  two  equations  yield 


or 


tan  e  = 


-2B  (Y-Cy) 

ic  (Z-GJ  +'  D 
z 


e  =  tan 


-1 


(-2B.  "'V _ ) 

\  2C  (Z-G  )  +  D  / 


(55) 


By  translating  and  rotating  the  plane  surface  from  the 
reference  plane  to  the  incident  point  of  an  input  ray  on  a 
quadric  surface,  the  plane  surface  system  can  be  used  to  trace 
a  ray  from  its  input  plane  to  the  output  plane  for  a  variety  of 
optical  element  shapes.  The  relationship  between  '-.he  plane 
surface  system  and  a  quadric  surface  is  summarized  as  follows: 


X1D  =  X10  +  (Z1D~RI) 

cos  a^/cos 

(56a) 

Y1D  *  Y1Q  +  (Z1D”V 

cos  6^/cos 

(56b) 
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,  </i +  0  and  B2  -  4^(? >  0 

-  £/  £  ,  e#  *  0  and  6  0 

n  -  sin"1  (2A  (X-Cx)/N) 


e 


/  “2B  (Y-Cy>  \ 
\  2C  (Z-G^)  +  D  / 


(56c) 

(56d) 

(56e) 
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Chapter  4 
RAYTRACE  PROGRAM 


DATA  VARIABLES 

The  ray  tracing  program,  RAYTRACE,  has  been  prepared  for 
running  on  the  CDC  3200  computer  using  a  Gould  electro-static 
printer  for  the  graphics  output.  Appendix  B  describes  and 
lists  each  routine  in  the  RAYTRACE  program  along  with  a  flow 
chart  of  the  program.  The  input  data  necessary  for  tracing  a 
beam  through  an  optical  surface  are  stored  in  the  DV  array. 

Each  time  a  beam  is  to  be  traced  through  an  optical  surface, 
the  DV  array  is  initialized  to  zero  except  for  the  run  number, 
the  magnitude  of  the  DV  array,  and  the  number  of  elements  in  the 
optical  system.  A  list  of  the  DV  array  variables  is  shown  in 


Most  of  the  variables  in  the  DV  array  are  self  explanatory. 

The  variables  from  DV{1)  to  DV(21)  are  the  parameters  required 
to  describe  the  optical  surface  for  the  simulation  of  a  beam 
passing  through  that  surface.  These  parameters  have  been 
described  in  the  previous  chapters.  DV(22)  describes  the  optical 
shape  of  the  surface.  If  the  variable  is  equal  to  one  the  shape 
is  convex.  For  a  concave  surface  the  variable  is  set  to  minus- 
one.  It  is  important  to  note  that  all  angular  input  data  must  be 
radians.  In  addition,  the  units  for  the  linear  data  must  agree. 

Only  one  input  parameter  DV(23)  is  needed  for  the  graphics 
output  of  a  cross  sectional  view  and  a  knife  edge  scan  at  the 
output  reference  plane  of  an  optical  surface.  This  parameter  is 
the  radius  of  the  cross  section  one  wishes  to  observe  at  the 
output  distance  from  the  surface  and  must  be  greater  than  the 
beam  radius  in  order  to  plot  the  data.  For  reasonable  scales 
of  the  axes,  this  radius  should  be  some  even  multiple  of  four. 

DV(24)  and  DV(25)  provide  boundary  conditions  for  the  beam  at  the 
output  reference  plane.  These  parameters  are  the  maximum  and 
minimum  radius,  respectively. 

PROGRAM  CONTROLS 

The  user  has  continuing  control  over  the  program  through  the 
use  of  the  console  typewriter  and  sense  switches.  After  a  beam 
has  passed  through  an  optical  surface,  the  run  number,  the  surface 
number,  and  the  number  of  ray  failures  are  printed  on  the  typewriter 
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Tabla  3 

DV  ARRAY  VARIABLES 


DV{1) 

DV(16) 

D 

DV(2) 

*X 

DV (17) 

E 

DV{3) 

*R 

DV(18) 

C 

X 

DV(4) 

X 

DV(19) 

C 

y 

DV<5) 

e 

DV(20) 

c 

z 

DV(6) 

n 

DV(21) 

F 

DV(7) 

X1D 

DV(22) 

Optical  Shape 

DV(8) 

yid 

DV(23) 

Radius  of  Cross  Sectioned. 

Area 

DV(9) 

Z1D 

DV(24) 

R 

max 

DV(10) 

RP 

DV(25) 

R  * 
lain 

DV(U) 

N 

DV (26) 

Number  of  Present  Optical 
Element 

DV(12) 

S 

DV{27) 

Run  Number 

DV(13) 

A 

DV{28) 

Magnitude  of  DV  Array 

DV (14) 

B 

DV (29) 

Number  of  Elements  in  Optical 
System 

DV(15) 

C 

DV(30)-DV (50) 

Blank 
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in  order  to  help  the  user  decide  on  what  course  of  action  to  take. 
If  any  graphics  output  is  desired,  sense  switch  (2)  is  set.  If 
any  graphs  are  plotted,  the  graph  can  be  replotted  immediately 
with  a  different  dimension  or  scale  by  an  appropriate  input 
through  the  typewriter.  Setting  sense  switch  (3)  will  list  the 
output  linear  and  angular  positions  of  each  ray  in  the  beam. 

The  beam  can  be  retraced  through  the  same  optical  system  by 
setting  sense  switch  (4)  and  inputting  the  new  parameters  through 
the  typewriter  for  the  appropriate  optical  surface.  Initially, 
the  input  beam  to  the  optical  system  has  been  generated  and  stored 
in  a  data  file  on  disc.  After  the  beam  has  been  traced  through 
the  system  it  may  be  desirable  to  transfer  this  beam  to  the  input 
beam  data  file.  This  can  be  accomplished  by  setting  sense  switch 
(5) .  As  one  becomes  familiar  with  the  program,  some  of  the 
messages  can  be  deleted  by  setting  sense  switch  (6) . 


DATA  STRUCTURE 


a. 


Input 


The  data  for  each  optical  system  are  read  from  the  card 
reader  and  stored  in  a  data  file.  The  data  for  each  optical 
system  consist  of  heading  information  and  data  for  each  successive 
surface  in  the  system.  As  the  beam  is  traced  through  a  surface 
the  data  for  that  surface  is  transferred  from  the  data  file  to 
the  DV  array.  Figure  11  shows  a  typical  setup  for  an  input  data 
deck.  If  the  same  input  beam  is  used  or  if  sense  switch  (5)  is 
used  to  transfer  a  resultant  output  beam  to  the  input  beam  data 
file,  the  input  data  deck  can  consist  of  more  than  one  optical 
system. 

The  heading  cards  for  an  optical  system  are  placed  in 
front  of  the  system's  surface  data.  This  information  is  read 
and  stored  in  memory  as  Hollerith  constants  from  column  one 
through  column  seventy-one.  As  each  card  is  read,  column  seventy- 
two  is  tested  for  a  zero,  a  one,  or  a  two.  A  zero  means  that 
more  heading  information  will  follow.  A  one  indicates  that  all 
heading  information  has  been  read  into  memory.  The  data  card  with 
a  one  in  column  seventy-two  must  separate  the  heading  cardo  of 
the  system  from  its  surface  data.  A  two  in  column  seventy-two 
signifies  the  end  of  the  computer  run  and  this  card  must  be  the 
last  card  of  the  input  data  deck. 

After  the  heading  information  has  been  read  as  indicated 
by  the  one  card,  the  data  of  the  first  optical  surface  of  the 
system  are  read  and  stored  in  the  DV  array.  The  format  of  this 
data  is  4(14,  E14.7)  for  each  card.  The  location  of  each  variable 
in  the  DV  array  is  denoted  by  the  14  format  and  the  value  of  that 
variable  by  the  E14.7  format.  The  magnitude  of  the  DV  array  is 
fifty.  As  each  data  card  is  read,  the  location  of  each  variable 
is  compared  with  fifty.  When  a  comparison  matches,  all  of  the 
input  data  for  that  particular  surface  have  been  read  and  stored 
in  memory.  Therefore, the  last  data  card  cf  the  input  data  for 
an  optical  surface  should  be  an  integer  fifty  in  14  format.  The 
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data  for  each  surface  of  an  optical  system  are  stacked  in  the  same 
order  as  they  occur  in  the  system.  If  a  bean  is  to  be  traced 
through  another  optical  system,  the  data  deck  is  set  up  as  outlined 
and  stacked  behind  the  previous  system  a3  shown  in  Figure  11. 

b.  Output 

The  first  piece  of  information  written  on  the  line  printer 
is  the  heading  information  for  an  optical  system.  Then  a  listing 
of  the  data  stored  in  the  DV  array  for  the  current  surface  is 
printed.  After  the  listing  of  the  DV  array,  two  rows  of  integer 
numbers  are  printed.  These  numbers  represent  the  number  of  rays 
stored  in  each  data  block  on  the  disc  before  and  after  the  rays 
have  passed  through  the  optical  surface.  Then  the  total  number 
of  rays  that  have  failed  to  pass  through  the  current  optical 
surface  is  listed.  If  sense  switch  (3)  is  set  the  output  beam 
data  for  the  current  optical  surface  are  listed  for  each  ray. 

The  format  is  (15,  1H) ,  4E20.7)  for  each  ray.  This  information 
is  the  coordinate  of  the  ray  in  the  output  plane  (X^D,  Y^D)  and 

the  angular  deviation  (j!.,  of  the  ray  from  the  output  reference 
ray.  *  * 

SAMPLE  PROBLEMS 

a.  Ellipsoidal  Mirror 


Several  problems  with  known  results  are  presented  to  test 
the  RAYTRACE  program  and  to  show  its  uses.  A  concave  ellipsoidal 
mirror,  as  shown  in  Figure  12,  is  chosen  as  the  first  sample 
problem  because  of  one  of  its  properties.  If  a  ray  originates  at 
one  focus  point  of  an  ellipsoidal  mirror  the  ray  will  reflect 
from  the  mirror  and  pass  through  the  other  focus  point.  The 
parameters  of  the  quadric  surface  are  a  *  6,  b  =  6,  and  c  ®  10. 

If  the  reference  system  of  the  quadric  surface  is  placed  at  the 
vertex,  the  center  of  symmetry  is  at  Z  =  10.  The  input  and  output 
reference  planes  are  at  the  focus  points  Z  *  2  and  Z  *  18, 
respectively.  The  input  reference  ray  is  coincident  with  the  Z-axis. 
With  this  sample  problem  in  mind,  a  beam  from  a  point  source  was 
generated  and  stored  in  the  input  beam  data  file.  With  the  point 
source  at  the  focus  point  Z  *  2,  the  beam  should  be  focused  at 
the  other  focus  point  after  tracing  the  beam  through  the  mirror. 
Figure  13  lists  the  input  data  deck  used  for  this  sample  problem. 

The  two  runs  have  the  output  reference  planes  at  two  different 
distances  from  the  surface.  Figure  14  shows  the  cross  sectional 
view  of  the  beam  focused  at  the  focus  point  of  Z  =  18.  When  the 
output  reference  plane  is  moved  to  Z  *  2,  Figure  15  shows  the 
beam  pattern  at  this  distance  from  the  mirror.  A  knife  edge  scan 
of  the  beam  at  the  distance  Z  *  2  is  shown  in  figure  16. 
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Figure  13  Input  Data  for  Ellipsoidal  Mirror 
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Figure  14  Cross  Section  of  Beam  at  Focus  Point 
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b.  Elliptic  Cylinder 

Since  the  RAYTRACE  program  is  not  limited  to  just  quadric 
surfaces,  an  elliptic  cylinder  can  be  formed  from  the  previous 
example  by  setting  either  a  =  0  or  b  =  0.  By  setting  a  =  0  the 
axis  of  the  cylinder  is  parallel  to  the  X-axis.  If  b  =  0,  then 
the  axis  of  the  cylinder  is  parallel  to  the  Y-axis.  Using  the 
same  parameters  of  the  previous  problem  the  focus  points  are 
Z  *  2  and  Z  =  18.  By  placing  a  point  source  at  the  focus  point 
Z  *  2  the  reflected  beam  should  be  focused  along  a  line  parallel 
to  the  axis  of  the  cylinder  at  Z  *  18  as  shown  in  Figure  17.  The 
input  data  are  listed  in  Figure  18  for  two  cylindrical  mirrors  of 
different  cylindrical  axes.  Figures  19  and  20  show  that  the  cross 
section  of  the  reflected  beam  in  the  plane  of  the  focus  point 
Z  =*  18  is  a  line  parallel  to  the  axis  of  the  cylindrical  mirror. 

c.  Telescope 

This  sample  problem  will  show  the  use  of  the  focusing 
parameter  or  DV(21).  This  problem  is  depicted  by  a  telescope 
which  consists  of  a  convex  spherical  mirror  and  a  concave 
ellipsoidal  mirror  as  shown  in  Figure  21.  The  reference  coordinate 
systems  of  the  two  surfaces  are  placed  at  the  vertices  and 
separated  by  a  distance  of  124.2.  The  spherical  mirror  is  surface 
one  with  the  center  of  symmetry  at  Z  =  -73.6.  Surface  two  is  the 
ellipsoidal  mirror  with  the  quadric  constants  a  =  b  =  675.13995 
and  c  **  1415.5712.  The  center  of  symmetry  is  at  Z  =  1415.5712. 

It  should  be  noted  that  the  dimensions  for  a  particular  surface 
are  measured  from  its  respective  reference  coordinate  system. 

A  source  of  parallel  rays  is  considered  to  be  somewhere  between 
the  two  mirrors.  With  the  separation  distance  of  124.2  the  beam 
will  focus  at  infinity  after  reflecting  off  the  two  surfaces. 

By  increasing  the  distance  between  the  mirrors,  the  beam  will 
converge  at  some  distance  less  than  infinity.  To  focus  the  beam 
at  the  far  field  distances  of  70,000  and  100,000  the  center  of 
symmetry  is  moved  by  -.396  and  -.263,  respectively.  This  in 
effect  increases  the  distance  between  the  mirrors.  The  input  data 
for  focusing  the  telescope  at  these  two  distances  are  listed  in 
Figure  22.  Figures  23  and  24  show  that  the  telescope  can  focus 
the  beam  at  the  far  field  distances  by  the  use  of  the  focusing 
parameter. 

d.  Telescope  with  Plane  Mirror 

Figure  25  shows  a  telescope  of  the  previous  problem  with 
the  addition  of  a  plane  mirror.  The  orientation  of  the  various 
coordinate  systems  needed  to  trace  a  beam  through  the  system  is 
shown  for  each  surface.  A  dot  represents  the  X-axis  out  of  the 
page  and  a  cross  represents  the  X-axis  into  the  page.  It  should 
be  pointed  out  that  the  output  of  surface  one  is  not  oriented 
with  the  input  of  surface  two  even  though  the  planes  of  incidence 
are  coincident.  Therefore,  a  rotation  of  tt  radians  is  necessary 
to  orient  the  two  coordinate  systems. 
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Figure  18  Input  Data  for  Cylindrical  Mirror 
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Figure  19  Reflected  Beam  from  Cylindrical  Mirror  Parallel  to  X-Axis 
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Figure  20  Reflected  Beam  from  Cylindrical  Mirror  Parallel  to  Y-Axis 
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Figure  22  Input  Data  for  Telescope 
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The  previous  example  has  shown  that  tha  beam  can  be  focused 
at  some  distance  in  the  far  field.  If  the  plane  mirror  has  some 
angular  deviation  from  its  reference  coordinate  system,  the  focused 
beam  in  the  far  field  will  move  about  the  origin  of  the  output 
reference  plane.  Figure  26  shows  a  listing  of  input  data  where 
e  »  n  ®  6  milliradians  for  the  plane  mirror.  Figure  27  shows  how 
far  the  focused  beam  has  translated  from  the- origin. 

e.  Two  Thick  Lenses 

So  far  the  examples  have  been  concerned  with  mirrors,  the 
RAYTRACE  program  can  handle  the  refraction  case  also.  Two  thick 
lenses  are  used  to  demonstrate  this  capability  as  shown  in 
Figure  28.  A  point  source  is  placed  at  the  focal  point  of  the 
first  lens.  After  the  beam  passes  the  first  lens  the  rays  should 
become  parallel.  On  passing  the  second  lens  the  rays  should 
focus  into  a  point  at  the  second  focal  point  of  the  system.  The 
listing  of  the  input  data  and  parameters  of  the  system  are  given 
in  Figure  29.  Figure  30  shows  the  beam  focused  at  the  second 
focal  point  of  the  system  after  it  has  been  traced  through  the 
two  lenses. 
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16 

00000000+00 

17 

00000000+00 

18 

00000080+00 

19 

00000000+00 

20 

00000000+00 

21 

00000000+00 

22 

00000000+00 

23 

16000000+01 

24 

16000000+01 

23 

29 

30 

1 

00008000+00 

30000000+00 

26 

10000000+00 

27 

70000000+00 

28 

30000000+01 

13707963+00 

2 

10016000+01 

3 

12000000+02 

4 

00000000+00 

3 

00000000+00 

6 

00000000+00 

7 

00000000+00 

8 

00000000+00 

9 

00000000+00 

10 

-10000000+00 

11 

10000000+00 

12 

10000000+00 

13 

10000000+00 

14 

10000000+00 

13 

10000000+00 

16 

00000000+00 

17 

3416'?  <.00+03 

18 

00000000+00 

19 

00000000+00 

20 

-73600000+01 

21 

-"7600000-01 

22 

10000000+00 

23 

40000000+01 

24 

40800000+01 

23 

30 

1 

60000000+00 

26 

20000000+00 

15707963+80 

2 

41990000+00 

3 

70000000+84 

4 

00000000+00 

3 

00000000+00 

6 

00000000+00 

7 

00000080+08 

8 

00000000+00 

9 

00000000+00 

10 

-10008000+00 

11 

10000600+00 

12 

-10000000+00 

13 

21938773-06 

14 

21938775-06 

15 

49904138-07 

16 

00000000+00 

17 

10000000+00 

18 

00000000+00 

19 

00000000+00 

20 

.1'  155712+03 

21 

00000000+00 

22 

-10000000+00 

23 

40000000+02 

24 

40000000+02 

23 

30 

00800000+00 

26 

30000008+00 

Figure  26  Input  Data  for  Telescope  with  Plane  Mirror 
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CROSS  SECTIONAL  VIEW 
RUN  NO. 7  5URFRCE  NO.  3 

DATE  12/18/72  »  D ISTRNCE  =  78000.00 


i 

Figure  27  Translation  of  Focused  Beam  Due 

Angular  Deviation  of  Plane  Mirror 
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RUN  NO.  08008  SURFACES  =*  i  THRJ  4  DATE  =  21  DEC  1973  0 
TWO  THICK  LENSES  0 
INPUT  BERN  IS  POINT  SOURCE  RT  FIRST  FOCAL  POINT  OF  SVSTEM  0 


OUTPUT  BERM  IS 

;  RT 

SECOND  FOCAL 

POINT 

OF  SVSTEM 

0 

A 

1 

13707963+80 

2 

33434343+00 

3 

10000000+00 

4 

X 

00000000+00 

5 

00000000+00 

6 

00000000+00 

7 

00000000+00 

8 

00000000+00 

9 

00800080+00 

10 

10000000+00 

11 

66666667-01 

12 

10006000+00 

13 

10000000+00 

14 

10000000+00 

15 

10000000+00 

16 

00800000+00 

17 

90000000+00 

18 

00000000+00 

19 

00000000+00 

28 

-30000000+00 

21 

00000000+00 

22 

10000000+00 

23 

40000000+00 

24 

40800080+00 

23 

29 

30 

1 

00000000+00 

40000000+00 

26 

10000000+00 

27 

80000000+00 

28 

30000000+01 

13707963+00 

2 

10000000+00 

3 

31818182+00 

4 

00000000+00 

5 

00000000+00 

6 

00000000+00 

7 

00000000+00 

8 

00000e00+00 

9 

00000000+00 

10 

10000000+00 

11 

15000000+00 

12 

-10000000+00 

13 

10000000+00 

14 

18000008+00 

13 

10000000+00 

16 

00000000+00 

17 

23008000+01 

18 

00000000+08 

19 

00080000+00 

20 

50000000+00 

21 

00000000+00 

22 

-10000000+00 

23 

48000000+00 

24 

40000000+00 

25 

30 

1 

00000000+00 

26 

28000000+00 

15707963+00 

2 

35454545+00 

3 

10000000+00 

4 

00000000+00 

5 

00000000+00 

6 

00000800+00 

7 

00000000+08 

8 

00000000+00 

9 

00000000+00 

10 

10000000+00 

11 

66666667-01 

12 

18000000+00 

13 

10000000+00 

14 

10000000+00 

13 

10000000+80 

16 

00000000+00 

17 

90000000+00 

18 

08000000+00 

19 

00000000+80 

20 

-30000000+00 

21 

00000000+00 

22 

10000000+00 

23 

40000000+00 

24 

40000000+00 

23 

30 

1 

00000000+00 

26 

30000000+00 

13707963+00 

2 

10000000+00 

3 

31818182+00 

4 

00000000+00 

5 

00000000+00 

6 

00000000+00 

7 

00000000+00 

8 

00000000+00 

9 

00000000+00 

10 

10000000+00 

11 

15000000+00 

12 

-10000000+00 

13 

10000000+00 

14 

10000000+00 

15 

10000000+00 

16 

00000000+00 

17 

25000000+01 

18 

00000000+00 

19 

00000000+00 

20 

50000000+00 

21 

00000000+00 

22 

-10000000+00 

23 

40000000+00 

24 

40000000+00 

23 

30 

00000000+00 

26 

40000008+00 
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Figure  29  Input  Data  for  Two  Thick  Lenses 
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Appendix  A 
GLOSSARY 


A  -  Parameter  in  the  basic  equation  of  a  quadric  surface 

a  -  Constant  used  in  quadric  equations 

-  Constant  formed  by  the  combination  of  terms  and  used  in  a 
quadratic  equation  for  solving  the  Z  coordinate  of  the 
incident  ray  on  a  quadric  surface 


B  -  Parameter  in  the  basic  equation  of  a  quadric  surface 

b  -  Constant  used  in  quadric  equations 

ft  -  Constant  formed  by  the  combination  of  terms  and  used  in  a 
quadratic  equation  for  solving  the  Z  coordinate  of  the 
incident  ray  on  a  quadric  surface 

C  -  Parameter  in  the  basic  equation  of  a  quadric  surface 

C  .  C  .  C  -  Center  of  symmetry  of  a  quadric  surface 
x  y  z 

c  -  Constant  used  in  quadric  equations 

C  -  Constant  formed  by  the  combination  of  terms  and  used  in  a 
quadratic  equation  for  solving  the  Z  coordinate  of  the 
incident  ray  on  a  quadric  surface 

D  -  Parameter  in  the  basic  equation  of  a  quadric  surface 

E  -  Parameter  in  the  basic  equation  of  a  quadric  surface 

P  ~  Parameter  which  focuses  a  quadric  surface  along  the  Z-axis 
of  symmetry 

G  ,  G ,  G  -  Parameters  formed  by  combination  of  terms  in  the 
x  ^  basic  equation  of  a  quadric  surface 

N  -  Ratio  of  nj/nR  for  a  refractive  surface  and  equal  to  one 
for  a  reflective  surface 

nj  •»  Index  of  refraction  of  the  medium  of  the  incident  ray 

nD  -  Index  of  refraction  of  the  medium  of  the  refracted  ray 

RP  -  A  constant  equal  to  one,  whose  sign  determines  the  direction 
of  the  output  ray  from  an  optical  surface  (+1  refractive 
surface,  -1  reflective  surface) 
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Rj  -  Distance  from  the  origin  of  the  input  reference  plane  to  the 
origin  of  the  plane  surface  reference 

R^  -  Distance  from  the  origin  of  the  plane  surface  reference  to 
the  origin  of  the  output  reference  plane 

S  -  A  constant  equal  to  one,  whose  sign  determines  the  concavity 
or  convexity  of  a  quadric  surface 

T-  -  Matrix  representing  the  transformation  of  vectors  from  the 
input  reference  plane  to  the  plane  surface 

T_  -  Matrix  representing  the  transformation  of  vectors  from  the 
plane  surface  to  the  output  reference  plane 

Vj  -  Incident  ray 

VR  -  Refracted  or  reflected  ray 

X. ,  Y . ,  z .  -A  right  handed  orthogonal  coordinate  system  where 
3  3  3  j  »  1,  2,  ...  n 

Xl'  Yl'  zi  “  Plane  surface  reference  coordinate  system 

X2»  Yj»  Z2  -  Input  reference  coordinate  system 

X3,  Y^»  Z3  -  Coordinate  system  attached  to  the  plane  surface 

X^,  Y^,  Z^  -  Output  reference  coordinate  system 

X._,  Y._,  "  Translation  of  a  coordinate  system  from  its 

3°  3  3  reference  position 

Xa0,  Y-0,  Z._  -  Coordinate  point  in  one  coordinate  system  in  terms 
jU  ju  Qf  another  coordinate  system 


X 


X 


Xj?'  YjPf  ZjP  “  Incident  point  of  ray  on  Z  =  0  plane 


8.,  y.  -  Angles  that  define  the  direction  of  a  ray  with 
3  3  respect  to  the  X,  Y,  and  Z  axes  and  their  cosines 

are  known  as  the  direction  cosines  of  the  ray 


e  -  Angle  of  rotation  of  the  plane  surface  about  the  X^  axis 

n  -  Angle  of  rotation  of  the  plane  surface  about  the  Yi  axis 

0  -  Angle  of  the  input  reference  ray  makes  with  the  negative  Y^ 

axis  in  the  plane  of  incidence 

0_  -  Angle  the  output  reference  ray  makes  with  the  positive  Y^ 

R  axis  in  the  plane  of  incidence 

\  -  Angle  between  the  positive  X,  axis  of  an  output  reference 

plane  and  the  positive  X0  axis  of  the  following  input  reference 
plane 
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-  Angle  of  incidence 
tr  -  Angle  of  refraction  or  reflection 

j g  -  Angular  deviation  of  the  input  ray  from  the  positive  Z2 
2  axis  in  the  *  0  plane 

0  -  Angular  deviation  of  the  output  ray  from  the  positive  z4 
4  axis  in  the  *  0  plane 

\| i  —  Angular  deviation  of  the  input  ray  from  the  positive  Z2 
2  axis  in  the  X2  =*  C  plane 

ij,  -  Angular  deviation  of  the  output  ray  from  the  positive  Z^ 
4  axis  in  the  X4  *  0  plane 


\.'\ .  •<  -'  '•'  ,....-.  ....  ,  ,,  ■'•  ■■'■■■  •  •■•■  ifi'i  ' -iAMf  -• 
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Appendix  B 

RAYTRACE  PROGRAM  DESCRIPTION 


A  brief  description  of  each  routine  in  the  RAYTRACE  program  is 
given  along  with  a  listing.  A  flow  diagram  of  the  program  is  also 
presented.  The  parameter  R  is  a  Hollerith  constant  containing 
the  date. 

RAYTRACE  -  The  main  routine  sets  up  the  input  beam  data  and 
the  surface  parameters  for  each  optical  system.  This  routine 
exercises  overall  control  of  the  program  options. 

BCOND  (RMAX2 ,  RMIN2,  NRY)  -  After  the  beam  has  passed  through 
the  optical  surface  to  the  output  reference  place,  each  ray  of  the 
beam  is  compared  with  a  set  of  boundary  conditions  for  a  circular 
beam  at  the  output  reference  plane.  RMAX2  is  the  square  of  the 
maximum  radius  and  RMIN2  is  the  square  of  the  minimum  radius.  NRY 
is  the  number  of  rays. 

BEAMIN  -  Beam  data  that  has  been  generated  and  stored  on  disc 
in  a  data  """’file  is  transferred  to  a  scratch  file.  The  number  of 
data  blocks  and  the  number  of  rays  in  each  block  are  stored  in  the 
NR  array. 

CROSSECT  (R)  -  The  data  necessary  for  plotting  a  cross 
sectional  view  of  a  beam  at  the  output  reference  plane  of  an 
optical  surface  is  set  up  and  stored. 

DATE  (R)  -  The  date  is  retrieved  from  a  system  subroutine  and 
is  stored  in  R  as  a  Hollerith  c  instant  in  the  form  MM/DD/YY. 

GRAPHS  (R)  -  This  subroutine  receives  the  generated  data  from 
CROSSECT  or  SCAN  and  plots  the  output  on  the  electrostatic  printer. 

INCIDENT  (NRY)  -  This  subroutine  calculates  the  incident 
position  of  the  input  ray  on  a  quadric  surface  and  the  angular 
rotation  of  the  plane  surface  system  required  to  place  the  plane 
tangent  to  the  incident  point.  NRY  represents  the  number  of  rays. 

INTRFACE  -  The  output  ray  parameters  are  interfaced  with  the 
input  ray  parameters  of  the  following  surface  when  the  output 
reference  system  is  out  of  phase  with  the  following  input  reference 
system. 

MATVEC  (A,  B,  C,  N)  -  This  subroutine  combines  six  different 
matrix-vector  or  matrix-matrix  operations  in  one  subroutine.  A, 

B,  and  C  are  3x3  matrix  arrays  stored  column  wise  or  vectors 
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depending  on  N  where  A,  B,  and  N  are  inputs  and  C  is  the  output. 

If  capital  letters  denote  3x3  matrices  and  lower  case  letters 
denote  vectors,  MATVEC  computes: 

N  =  0,  c  *  Ab 

N  »  1,  C  =  ATb 

N  =  2,  C  =  AB 

N  -  3,  C  =  ATE 

N  *  4,  C  =  ABT 
T  T 

N  =  5,  C  =  AXB 

where  T  represents  the  transpose  of  a  matrix. 

READ  -  Prior  to  tracing  a  beam  through  an  optical  system,  the 
surfaces  are  read  and  stored.  An  entry  through  LIST  will  print 
the  parameters  of  the  particular  surface. 

RETRACE  -  If  it  is  desired  to  change  certain  pa.-ameters  in  an 
optical  system  and  retrace  the  beam  through  the  system,  this 
subroutine  will  accept  input  from  the  typewriter  for  changes  to 
each  surface  and  will  reload  the  input  beam. 

SCAN  (R)  -  The  data  necessary  for  plotting  a  knife  edge  scan 
of  a  Hearn  at  the  output  reference  plane  of  a  surface  is  set  up  and 
stored. 

SINCOS  1  (A,B)  -  The  sine  and  cosine  functions  of  angle  A  in 
radians  are  calculated  with  the  output  stored  in  the  B  array. 

B (1 )  contains  the  sine  of  A  with  the  cosine  of  A  in  B(2). 

TlxACE  (NRY,  AZERO)  -  A  number  of  rays,  NRY,  are  traced  from 
the  input  reference  plane  through  the  plane  surface  system  to  the 
output  reference  plane.  If  =  tt/2,  then  AZERO  equals  0.  For  all 

other  values  of  AZERO  equals  1. 

TRANSFER  -  If  it  is  desired  to  save  the  output  beam  data  for 
later  use,  this  subroutine  will  transfer  the  beam  data  from  a 
scratch  file  to  a  saved  data  file. 

WRITE  -  This  subroutine  will  list  the  output  beam  parameters 
on  the  line  printer. 
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Figure  B-l  RAY TRACE  Flow  Diagram 
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Figure  B-l  RAYTBACE  Flow  Diagram  (Cont.) 
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Figure  B-l  RAYTRACE  Flow  Diagram  (Cont. ) 
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c 

C 

c 

c 

c 

c 

c 

c 
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c 

c 

c 

c 

c 

c 

c 
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c 
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c 

c 

c 

c 

c 
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PROGRAM  KRYlkMlUL 

RAYTRACE  IS  THE  Mil  IN  ROUTINE  I  UR  IRflClNU  RAYS  II  (ROUGH 

AN  OPTICAL  SYSTEM. 

9  JANUARY  19?4  <RKVISFD> 

DATA  PILES  ON  DISC  : 

BEAM  THIS  Pli  K  CnN  LA  i  NS  A  GENERA  I  PD  BhiilT  WHERE  EACH  RAY  IS 

DPP  I  NED  BY  (WO  I  I  NEAR  POSIT  IANS  <X,  Y>  AND  TWO  RNQUL  flR 
POSITIONS  < PH I;  PSD.  (HIS  I  IIP  HA".  A  MAXIMUM  OF  3  DATA 
Bl  OCKS  WITH  A  MAXIMUM  UP  200  RAYS  IN  F.ACH  BLOCK. 

NRYF  —  THIS  FILE  HAS  1  HE  NUMHPR  OP  DlliA  BLOCKS  STORED  IN  THE 

BEAM  PUP  AND  THE  NUN  APR  OP  kllYS  i  N  EACH  BLOCK. 

RHYF  SCR  A I  Cl  I  PUP:  USI  D  P  UR  RAY  DATA  AS  I  HP  BEAM  IS  TRACED 
THROUAH  I  HE  SYS  I  EM. 

XQRF  •-  TUP.  X  POSITION  POR  A  PllRTILUIllR  GRAPH  IS  STOkED 
IN  THIS  PILE. 

YQRF  --  THE  Y  POSITION  IOR  A  FliRTloiJI.  MR  OklTPH  IS  STORED 
IN  THIS  FII.P. 

STOR  —  THE  PAkAMPT  PRS  up  EACH  OPTlCHl  SURFACE  ARE  STOkED 
IN  THIS  PIl.P. 

ARRAYS  IN  THE  RAYTkACE  PROmRAM  : 

RAY  --  AS  INPUT  THIS  ARRAY  CONTAINS  I  HE  I  I NEAR  AND  ANGULAR 
POSITIONS  OP  EACH  RaY  PROM  A  I'ATII  BLOCK  AT-  THE  INPUT 
REFERENCE  HANP.  AS  01)  1  PUT  HUS  ARRAY  CONTAINS  THE  RAY 
POSITIONS  A)  THE  OUlPUf  REP PRINCE  PLANE. 

DEV  —  THIS  ARRAY  CoNIAINS  HIP  ANAUl  Ilk  <  I  PS  I  LON;  ETA)  AND 

t. INPllR  <X1D;  YiD;  XI D)  DP.VlmlnAS  01  A  PI  ONE  SURPACE  PROM 
THE  SURE  ACE  REP  I  REACH  PLANE  To  Till  INCIDENT  POINT  OF  AN 
INPUT  RttV  ON  A  DU  ADR  If:  SUKI  HOP 

TRF  THIS  ARRAY  1DP7N I  IF  IPS  RAY  Fit  1 1  ARPS  BY  SETTING  IRF  ~  1 

FoR  A  FAILED  RAY. 

DV  ---  THE  PARAMETERS  OF  A  PlIKflUMnR  OPlICAL  SURFACE 
ARE  CONTAINED  IN  THIS  ARRAY. 

NR  --  THE  NUMBER  UP  RAYS  IN  EACH  Dll  I H  Bl  OCK  ARE  STORED 
IN  THIS  ARRilY. 

PROGRAM  OPTIONS  : 

SET  SENSE  SWITCH  <?)  -  GRAPHICS 

SET  SENSE  SWITCH  <3)  -  I  1ST  UU 1 PU f  HHiM  PnRAMETERS 

SET  SENSE  SWITCH  <4>  -  RETRACE  BEAM  THRU  SYSTEM  WITH  NEW  INPUT 

SET  SENSE  SWITCH  <5>  -  TRANSFER  BP AM  PARAMETERS  TO  BEAM  FILE 

SET  SENSE  SWITCH  <6)  -  DHL!  fh  MOS I  IViEWRITTFN  MESSAGES 

COMMON  RAY <4, 250); DEV<b, 256),  IRF <250);  DV<50>; NBI  K; NR<24> 

2  FORMAT  <  1H1  ) 

4  FORMAT  <  71H 

1  ID 


s 
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c 

c 

c 


c 

c 

c 


6  FORMAT  <  1H@, 41  I  IDO  NOT  SET  SENSE  bWIiCiFs  UNlIL  REQUESTED  ) 

8  FORMAT  <  1018  > 

10  FORMAT  <  //VSX23HNUI1BER  OF  RAV  FAILURES  -  ,14//  > 

12  FORMAT  <  1H8, 13HRUN  NUMBER  =  ,  15/ 

125H0PTICAL  SURFACE  NUMBER  =  ,  14/ 

225HNUMBER  OF  RAV  FAILURES  =  , 14  > 

14  FORMAT  <  1H@,  L8H0RAPHICS  C  /  20HSF.T  SENSE  SWITCH  <2>  1 

15  FORMAT  <  1H0,  42HSET  REMOTE  SWITCH  ON  FLFC I ROS  TATI  C  PRINTER  / 
132HTVPE  10  FOR  CROSS  SECTIONAL  VIHWi'Si'Hl VPE  01  FOR  KNIFE  EDGE  SCAN 
2/23HTVPE  11  FOR  BOTH  GRAPHS  > 

15  FORMAT  <  1H0, 24HIVPE  IN  GRAPH  T  OS  CHOICE  > 

1 ?  FORMAT  <  211  > 

18  FORMAT  <  1H@,  2VHI  1ST  GUTPU 1  BEAM  PARAMF.'l  I  RS,  1HC  / 

1  20HSET  SENSE  SWIICH  <31  > 

20  FORMAT  <  1H@, 3?HTRANSFER  BEAM  PARllME  I ERS  TO  BEAM  FILE  ,  1HC/ 

1  20HSET  SENSE  SWITCH  05)  > 

2?  FORMAT  <  1H0, 41HKE TRACE  BEAM  THRU  SVsihM  WI1II  NEW  INPUT  C  / 

1  20HSET  SENmF  SWITCH  <41  1 


READ  AND  NRIIE  HEAD I NO 


CALL  DATE  <R> 

28  PRINT  2 
30  READ  4, I 
PRINT  4,  I 

IF  <  I  -•  1  >  30,  32,  60 


READ  IN  BEAM  AND  THE  OPIICAL  SURFACE  PARniiEl  ERS. 


32  CALL  BEflMIN 
CALL  READ 

33  IK  =  0 


READ  IN  PRESENT  OPUCAI  SURFACE  PARAMETER'-.. 


3<1  IK  =  IK  +  1 

CALL  FRF  AD  <  4HS  I  OR,  IK,  DV,  :1  00  > 

92  IF  <  I RS  I'ATF <4HST0R >  .  NF.  1  )  GO  10  99 


C 

L 

C 


PRINT  PRESENT  OPIIClH.  SURFACE  I'ARflMi-.  l  ERS.. 


CALL  LIST 
A2.ER0  --  0.  0 

IF  <  DV(1)  NS.  n.D/’0?963  1  A/ERO  -  1.(1 
IF  <  SSHTCHI  <61  .  KQ.  2  >  WRITE  <09,6  > 
PAUSE 

PR’NT  8.  <  NR<K),  K-  1,NBI  K  > 

N’F  =  0 


TRACE  EACW  bl  OCK  OF  WAVS  IMROUull  IMF  uFTIi.AL  SURFACE 
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DO  48  K  =  1*  NBLK 
NRY  =  NR<K> 

IF  <  NR<K>  .  EQ.  0  >  00  TO  48 
DO  36  I  -  1, NRY 
36  IRF< I >  =  0 

CALL  FREAD  <  4HRAYP, K; RAY, 8*NRY  > 

38  IF  C  IRS  fATF<4HRAYF)  .  NE.  1  )  00  10  38 

IF  <  DV<22>  .  NF.  0.  0  >  CALL  INCIDENT  <  NKV  > 

COLL  TRACE  <  NRY, A2ER0  ) 

CHECK  BOUNDARY  CONDI i IONS. 

RMAX2  =  DV<24)*DV<24) 

RMIN2  DV<25)*DV<25) 

CALL  6C0ND  <  RMAX2, RMIN2, NRY  > 

TOTAL  RAY  FA I I URFS. 

N  =  0 

DO  44  I  =  1, NRY 
IF  <  IRF  < I )  .  F.O.  0  )  00  TO  40 
HRF  =  NRF  +  1 
00  TO  44 
40  N  =  N  +  1 

DO  42  J  s  1,4 
4?  RAY< J, N)  -  RAY< J. I) 

44  CONTINUE 
NR<K)  «  N 

IF  <  NR<K>  . FQ.  0  >  00  10  48 

RESTORE  RAY  PARmMET  hRS  T  HAT  HAVE  MET  hOUNDARY  CONDITIONS. 

CALL  FWRITE  <  4HRAYF. K, RAY.  8*N  ) 

46  IF  <  IMS  I ATF <  4HRAYF )  .  NF.  X  )  00  TO  46 
48  CONTINUE 

PRINT  8,  <  NR<K>, K“l, NBLK  ) 

PRINT  10.  NRF 
NRUN  =  DV<2?) 

NPDV  =  DV<26) 

WRITE  <59,12)  NRUN, NPDV. NRF 
PROGRAM  01' 1  IONS. 

IF  <  SSW I  CHF<6)  .  F.Q.  2  >  WRITE  <59,14) 

PAUSE  2 

IF  <  SSWTCHF' < 2 )  .  FM.  2  )  00  TO  49 
IF  <  SSW!CHF<6>  .  FQ.  2  )  WRITE  <59,15) 

WRITE  <59,16) 

READ  <58, 1?)  13, 12 
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IF  <  II  .  EQ.  i  >  CflLI.  CROSSECT  <  R  > 

IF  <  12  .  EQ.  X  >  CALL  SCAN  <  R  > 

49  IF  <  SSW'i  CHF<6)  .  EQ.  2  )  WRITE  <59-18) 
PAUSE  3 

IF  <  SSWTCHF<3>  .  EQ.  1  >  CALL  WRITE 
IF  <  DV<4>  .  NE.  8.  0  >  CALL  INIRFmCE 
IF  <  DV<26)  ~  DV<29)  )  34-50-60 

50  IF  <  SSWYCHF<6)  .  EQ.  2  )  WRITE  <59-22) 
PAUSE  4 

IF  <  SSWTCHF<4>  .  EQ.  2  )  GO  TO  52 
CALL  RETRACE 
00  TO  33 

52  IF  <  SSWTCHK<6)  .  EQ.  2  >  WftI TE<59-  28) 
PAUSE  5 

IF  <  SSW'i CHF<5>  .  EQ.  1  >  CALL  TRANSFER 
00  TO  28 
68  CONTINUE 
END 
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SUBROUTINE  BCOND  C  RMRX2>  Kt1IN2<  NRV  J 
C 

C  31  OCTOBER  19?3 

C 

C  THIS  SUBROUTINE  WILL  TEST  EHCH  RRY  RURINSf  R  SET  OF  BOUNDARY 

C  COND11 IONS  <RMRX> RMIN>  FOR  R  CIRCUIRR  BERN  RT  THE 

C  OUTPUT  REFERENCE  PLANE. 

C 

COMMON  RRY<4, 250>,  DEV<5, 258>,  IRFC250), DV<S0>, NBLK, NR<24> 

DO  40  J  =  1, NRY 

IF  <  I RF  <  J )  .  EO.  1  >  00  TO  40 

X2  =  RRY<1j  J>*RRV<1> J) 

Y2  =  RRY<2, J>*RAY<2, J> 

R2  a  X2  +  Y2 

IF  <  R2  .  OT.  RMRX2  .OR.  R2  .  LT.  RMIN2  >  IRFCJ)  =  1 
40  CONTINUE 
RETURN 
END 
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SUBROUTINE  BERN IN 
C 

C  29  OCTOBER  1973 

C 

C  THIS  SUBROUTINE  WILL  READ  IN  BERM  DATA  THAT  HRS  RLRERDV  BEEN 

C  GENERATED  AND  STORED  ON  DISC  IN  THE  BERM  FILE.  EACH  RRV  OF  THE 

C  BERM  IS  DESCRIBED  BV  TWO  LINEAR  POSITIONS  <X,Y>  RND  TWO  ANGULAR 

C  POSITIONS  (PHI*  PSD.  THERE  IS  A  MAXIMUM  OK  250  RAYS  PER  BLOCK. 

C 

C  IN  ORDER  THAT  THE  GENERA TED  BEAM  IS  No T  DESTROYED  THE  BEAM  DATA 

C  IS  TRANSPEREO  FROM  THE  BEAM  FILE  TO  THE  RftYF  FILE. 

C 

C  THE  NRYF  FILE  CONTAINS  THE  NUMBER  OF  BLOCKS  STORED  IN  THE  BEAM 

C  FILE  AND  THE  NUMBER  OF  LAYS  IN  EACH  BLOCK. 

C 

COMMON  RAY <4/ 250), DEV<5, 250), IRF<250>, DV<50>, NBLK, NRC24) 

CALL  FREAD  <  4HNRYF, 1, NBLK, 25  ) 

30  IF  <  IRSTATF<4HNRYF)  .  NE.  1  )  GO  TO  20 
DO  40  K  =  1,  NBl  K 
NRY  =  NR<K> 

CALL  FREAD  <  4HBEAM,  K,  RAY,  S*NRY  >.  ’ 

32  IF  <  I RSI ATF C4HBEAM)  .  NE.  1  )  GO  TO  22 
CALL  FWRITE  <  4HRAYF, K, RAY, 8*NkY  > 

34  IF  <  I WS  T ATF < 4HRAYF )  .  NE.  1  )  00  TO  24 
40  CONTINUE 
RETURN 
END 
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SUBROUTINE  CKOSSF.CT  C  R  ) 


C 

C- 

C 

C- 

C- 

C- 

C 


10  JANUARY  1974  (REVISED) 


---  THIS  SUBROUTINE  STORES  ON  DISC  IN  THE  XGRF  FILE  ftND  THE  YQRF  FILE 
—  THE  DATA  NECF.SSARV  TO  PLOT  A  CROSS  SECTIONAL  VIEW  OF  A  BEAM 
— -  AT  THE  OUTPUT  REFERENCE  PI. ONE  OF  ANY  PARTICULAR  SURFACE. 


COMMON  RAY<4, 250), DEV<5, 250),  IRF(250>,  DV<50>,  NBLK,  NR<24) 
COMMONZPLTS/VIBUF  <  6000),  A(8),  P(3>,  X(8),  Y(8),  l.NE<2),  LI  NX/  LABL<5,  8) 
1  , XP<2000>,  XMRX,  XMIN,  DX,  YMAX,  YMIN,  DY 

DIMENSION  CC250, 2) 

EQUIVALENCE  (DEV, C) 

FORMAT  (20HCROSS  SECTIONAL.  VIEW) 

FORMAT  ( 20H  Y  AXIS) 

FORMAT  (  1H®, 30HTYPE  1  FOR  YES  ,  TYPE  0  FOR  NO  / 

1  40HREPLOT  CROSS  SECTION  WITH  NEW  DIMENSION!'  ) 

8  FORMAT  (  11  ) 

10  FORMAT  (  1H9/26HIF  MISTAKE  ON  INPUT  TYPE  1  7 
1  42HNEW  DIMENSION  X  :XXXXXXXX:XX  ,  Mr STAKE  X  > 

12  FORMAT  (  E12.  7,  11  ) 

14  FORMAT  (  1H@,  26HINPUT  M  IS  I AKE  RFITYPE  INPUT  ) 


2 

4 

6 


IX, 1H0  ) 


C 

C- 

C 


INITIALIZE  THE  GRAPH  PARAMF.  IFRS  FOR  THE  CROSS  SECTIONAL  PLOT. 


ENCODE 

ENCODE 

P(l) 

Y(7) 

LNE(l) 


XMAX  =  YMAX  = 


(  20,  2,  L  ABL(1,  1) 

(  20,  4,  l.ABI.d,  8) 

-  5.  00  t  P(2) 

*  Y(8) 

*  LNE(2) 
DV(23) 


«  4.  50 
=-2 


8.  00 
6.  25 
3 


XMIN  =  YMIN  =  -DV<23) 


DX 


=  DY 


DV(23)/4. 0 


C 

C- 

C- 

C 


— -  FROM  EACH  DATA  BLOCK  IN  RAYF  FILE  STORE  I  HE  I  I NEAR  POSITIONS 
—  -  OF  EACH  RAY  IN  XORF  AND  YGRF  FILE. 


DO  30  J  -  1, NBLK 
NRY  =  NR< J) 

IF  (  NRY  .  EQ.  0  )  GO  TO  30 
CALL  FRFAD  (  4HRAYF,  J,  RAY,  8*NRY  > 

22  IF  (  IRS  I  RTF (4HRAYF >  .  NE.  1  )  GO  TO  22 
DO  24  I  «  1, NRY 
DO  24  K  =  1,  2 
24  C( I, K)  =  RAY(K, I > 

CALL  FHR1TE  (  4HXGRF  ,  J,  C ( 1 , 1 ),  2rNRY  ) 
26  IF  (  I WS I  RTF' ( 4HXQRF  >  .  NE.  1  )  GO  10  26 
CALL  FHPITE  (  4HYGRF  ,  J,  C( 1 , 2 ), 2  +  NRY  ) 
28  IF  (  IWS  t  RTF (4HYGRF)  NE.  1  )  00  TO  28 
30  CONTINUE 
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32  CALL  GRAPHS  <  R  > 

—  OPTION  TO  REPLOT  GRAPH  WITH  NEW  SCALE. 

WRITE  <59,6  ) 

Cc  AD  <58, 8)  I 

IF  <  I  .  EQ.  0  >  GO  TO  38 

WRITE  <59,10) 

34  READ  <58,12)  DIM, MISTAKE 

IF  <  MISTAKE  .  EQ.  8  )  GO  TO  36 
WRITE  <59,14) 

GO  TO  34 

36  KM AX  =  VMAX  -  DIM 
XMIN  =  VMIN  -  -DIM 
DX  =  DV  =  DIM/4.  0 
GO  TO  32 
38  RETURN 
END 
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SUBROUTINE  DflTECR J 
C 

C -  THIS  SUBROUTINE  Cfll  l.S  THE  DATE  FROM  ACCOUNTS  AND  STORES  IT  IN  R 

C - AS  A  HOLtERITH  CONSTANT  OF  THE  FORM  MM/DD/'VY. 

C 


MACRO  <a> 

EXT 

ACCOUNTS 

LDAQ 

ACCOUNTS 

SHQ 

-6 

SHAQ 

-6 

SRCH 

X+3 

SHA 

-s 

SHR 

12 

ADAQ 

X 

UJP 

*■*•2 

BCD.  C 

ENDM 

REAL  A 

R  =  R 

RETURN 

END 

8.  08/00/'t30 
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SUBROUTINE  GRAPHS  C  R  5 
C 

c -  24  NOVEMBER  1973  <REVISED> 

C 
C 
C 
C 


-  THIS  SUBROUTINE  WIt.L  PLOT  DATR  GF.NERA1  ED  OR  STORED  BV  CROSSECT 

-  OR  SCAN  SUBROUTINES. 


4H 

,  4H 

4H 

4H 

,  4H 

4H 

4HNO. 

,  4H 

4H 

4H  = 

,  4H 

4H 

4H 

,  4H 

4H 

4H 

.  4H 

X  , 

4HAXIS)  , 

3  5- 

3<0. 

30), 

2<5.  O0>,  i. 

50,  5.  73,  4.  75 

6  )  * 

9.  73 

,  2<9 

».  50,  9.  25), 

9.  25  )  , 

S  )  = 

7<0. 

00), 

90.  0  )  ,  < 

C  LINX  -  4H 

C 

C- 

C- 

C 


COMMON  RAY<4, 250), DEV<5, 250), IRF<250), DVC50), NBLK, NR<24> 
COMMON/PLTS// IBUF <6060),  A<8>,  P<2),  X<8>,  Y<8),  I  NE<2>, LINX.  LR6L<5.  8) 
1  , XP<2000>» XMAX,  XMIN,  DX,  YMAX,  VMIN,  DV 

DIMENSION  YP<2000> 

EQUIVftl.F.NCE  <RAY, YP)  ,  <l.ABL<26),  RR) 

DHTfi/PLTSX<  <l.ABL<I),I  6,25  )  = 

14HRUN  , 4HN0. 

14HDATE, 4H 
14HSURF,  4HACE 


14H  ,  4H 

1  <  <  X<  I  >,  I 
1  <  <  V <  I  ),  I 
1  <  <  R<  I  ),  I 
2  FORMAT  <  4H2  —  0 ) 
RR  =  R 
NP  =  1 


—  -  READ  IN  1  HE  l. INEAR  POSITIONS  FOR  A  PAk  riLULAR  GRAPH 

—  AND  STORE  IN  THE  XP  AND  YP  ARRAYS. 

DO  26  J  =  1,  NBLK 
NRY  =  NR<  J ) 

IF  <  NRY  .  EQ.  0  )  GO  TO  26 

CALL  FREAD  <  4HX0RF, J, XP<NP), 2fNRY  ) 

22  IF<  IRSTATF (4HXGRF)  .  N!:..  1  )  GO  TO  22 
CALL  FREAD  <  4HYGRF, J, YP< NP ) . 2* NRY  > 

24  IF  <  IRS  f  ATF <4HYGRP )  .  NE.  1  )  GO  TO  24 
NP  =  NRY  +  NP 


IF 

<  Np-*-250  .  GE 

.  1998 

) 

GO  TO 

27 

26 

CONTINUE 

27 

NP 

=  NP  -  1 

IF 

<  NP  .  EQ.  8 

)  GO  TO 

32 

DO 

28  I  =  1, NP 

IF 

<  XP<I)  .  LT. 

XMIN 

) 

XP<I) 

=  XMIN 

IF 

<  XP<I>  .  GT. 

XMAX 

) 

XP(I) 

=  XMAX 

IF 

<  YP  Cl)  .IT. 

VMIN 

) 

YP  <  I  ) 

=  YU  IN 

IF 

<  YP  < I )  .  GT 

YMAX 

> 

YP<  I  ) 

=  YMAX 

23 

CONTINUE 

..  - 

PLOT  THE  GRAPH. 
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CALL  GRAPH I NT  <  IBUF<i>/ IBUF<6000>  > 

K1  *=  NP  +  1  *  K2  =  NP  +  2 

XP<K1>  =  XMIN  $  XP<K2>  =  DX 

YP<K1>  =  YMIN  *  YP<Kf:>  =  DY 

CALL  AXISE  (  0.  50/P<i>/LINX/  -4,  8.  00,  0.  0,  XP<Ki>,  XP<K2>  > 
CALL  AXISE  <  4.  50, 1.  90,  LINX.  4,  P<2>,  90.  ,  YP<Ki>,  YP<K2>  > 
XP<Ki>  =  XP<K1>  -  0.  30*XP<K2> 

YP<K1>  =  YP<K1>  -  1.  00*YP<K2> 

CALL  LINEE  <  XP/ YP,  NP/ 1/ LNF.<1>,  LNE<2>  > 

DO  30  J  =  1,8 

30  CALL  SYMBOLE  <  X< J>,  Y<  J)/ .  14/  l.ABL<i,  J>/ A<  J),  20  > 

CALL  NUMBERE  <  1.  50/  9.  50/  .  14/  DV<2?>z  0.  0/  -i  > 

CALL  NUMBERE  <  6.  70/ 9.  50/ .  14,  DV<26>/ 0.  0/ -1  ) 

CALL  NUMBERE  <  6.  50/ 9.  25/ .  14,  DV<3>  ,0.0/  2  > 

WRITE  <11/ 2) 

CALL  GRAPHOUr  <  11  ) 

32  RETURN 
END 
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SUBROUTINE  INCIDENT  <  NRY  ) 

C 

-  24  NOVEMBER  ±9?3  <REVISKD> 

-  GIVEN  AN  INPUT  RAY  DEFINED  BV  LINEAR  AND  ANGULAR  POSITION 

-  WITH  RESPECT  TO  AN  INPUT  REFERENCE  PLANE.  THIS  SUBROUTINE  WILL 

-  CALCULATE  THE  INCIDENT  POSITION  OF  THE  INPUT  RAV  ON  ft  QUADRIC 

C -  SURFACE  AND  THE  ANGULAR  ROTATION  OF  THE  SURFftCE  REFERENCE  PLANE 

C -  REQUIRED  TO  PLACE  THE  PLANE  TANGENT  TO  THE  INCIDENT  POINT. 

C 

COMMON  RAY<4/230>.  DEVO,  250),  IRF<250>#  DV<50>/  NBLK/  NR<24> 

DO  35  J  =  1, NRY 

TANPHI  =  TANF  <  RAY<3iJ>  i 

TANPSI  =  TANF  <  RAY<4, J>  > 

TPHI  =  -TANPHI  $  TPS I  =  TANPSI 

DELX  =  RAY(1<  J>  *  DEl.Y  =  -RAS'<2,  J> 

XLAM  =  DELX- DV<18> 

YLAh  =  DELY- DV<19> 

2LAM  =  DV<20)+DV<21) 

R1  =  DV<13)*TPHI  *  R2  =  DV<i4>*TPSI 

R3  *  Xl.AM-DV<2>*TPHI  *  R4  =  YLAM-DV<2>rTPSI 

R5  =  DV<i5WLAM 
C 

C - SOLVE  FOR  ?.  INCIDENT  POSITION. 

C 

ft  =  R1*TPHI  +  R2*TPSI  +  DV<15> 

B  =  -<  2.  0*<  R1-SR3  ♦  R2*R4  -  RS  >  +  DV(16>  > 

C  =  DV<13>*R3*R2  +  DV<14>+R4*R4  *  2LAM*<  R5  DVO!6>  >  -  DV<i?> 

Q  =  B*B  -  4.  G^-A-kC 

IF  <  Q  .  IT.  0.  0  )  GO  TO  38 

IF  <  A  .  EQ.  0.  0  .  AND.  B  .  EQ.  0.  8  >  GO  TO  38 

IF  <  A  .  Eft.  0.  0  ■>  GO  TO  28 

DEV<5»  J)  =  <  B  +  DV<12>>rS0RTF  <  Q  )  >/<  2.  I  *ft  ) 

GO  TO  25 

20  DEV<5<  J >  =  C/B 

25  R1  a  DEV<5, J)  -  0V<2> 

C 

C -  CALCULATE  X  AND  Y  INCIDENT  POSITIONS. 

C 

DEV<3#  J)  =  XLAM  +  R1*TPHI 

DEVC4.J)  =  YLAM  +  Rt^TPSI 

C 

C -  SOLVE  FOR  ROTATION  ANGIES  EPSILON  AND  ETA. 

C 

R1  =  2.  0*-DV<J  3>*<  DEV< 3/  J >  -  DVOG>  > 

R2  =  2.  0tDVO4>*<  DEV<4»  J)  -  D V09>  > 

R3  =  ?.  0i'DV<15)t'<  DE  V<5»  J )  -  21AM  )  +  DV<16> 

R4  =  SQP.TF  <  RlfRl  +  R2*R2  ♦  R3*kX  > 

DEV < 1»  J )  =  ATANK  <  -R22R3  > 

DEV<2»  J)  =  AS  INF  <  DV(22)fR12R4  ) 
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GO  TO  35 
38  IRF< J>  «  1 
35  CONTINUE 
RETURN 
END 
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SUBROUTINE  INTRFACE 
-  24  NOVEMBER  1973  <REVISED) 

-  WHEN  THE  INCIDENT  PLANES  OF  TWO  OPTICAL  SURFACES  ARE  NOT  PARALLEL 

- I  HIS  SUBROUTINE  INTERFACES  THE  OUTPUT  RAY  PARAMETERS  WITH  THE 

-  INPUT  RAY  PARAMETERS  OF  THE  FOLLOWING  OPTICAL  SURFACE. 

COMMON  RAY <4/ ?50>, D£V<5, 250>, IRF<250>,  DV<58>, NBLK, NR<24> 

DIMENSION  A<4)/  B<2> 

CALL  SINCOS1  <  DV<4>,B  ) 

DO  28  K  =  1. NBLK 

n  =  n:uk> 

IF  <  N  .  EQ.  0  >  Uo  TO  28 
CALL  FRSAD  <  4HRAYF, K, RAY. 8*N  ) 

18  IF  <  IRSTATF <4HRAYF>  .  NE.  1  >  GO  TO  18 
DO  24  J  =  1,  N 
DO  20  I  =  1,  2,  2 

A(  I  >  =  RAY<  1>  +  RAY<I+i, J>+B<1> 

20  ACI+1)  -  -RAY<  h J)*B(1)  +  RAYCI+i, J>*B<2> 

DO  22  I  =  1,  4 
22  RAY< I»  J)  =  A <  I  > 

24  CONTINUE 

CALL  FWRITE  <  ‘H1RAYF,  K,  RAY,  8*N  > 

26  IF  <  IWSi ATF<4HRAYF>  .  NE.  1  >  GO  TO  26 
28  CONTINUE 
RETURN 
END 
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biitjRoin  im-.  iiiii  v:  i;ci!,  it,  i  .  ro 

DIMENSION  HOO,  CCS),  i  <9>,  |kv^ 

IF  <N>  10,6,10 

10  Ou  10  <i>,  6,  li,  6.  0>,  N 
6  DO  61  J- 1, 3 
01. 

00  TO  70 
‘i  M2-1 

DO  36  K  1, 3 
Kt«Ki C 

DO  36  J-  K# K1 ,  s 
I-CM2)  tKd) 
ui  llx=M7!l 

70  II-  CM  1  )  Vi,  71  ■  vx 
VI  M1"l 

DO  ?.i  J  1 , 3 
VI  0(J>-0<J> 

00  TO  00 
V7.  I'M”? 

00  TO  <  VI,  ?<\,  V‘t,  VO,  VO),  N 
VI  DO  76  J-1 , 9 
V6  0< J>"H< J) 

00  10  00 
Vi  112=1 

DO  66  K-  1,3 
K1=Ki6 

DO  66  J'  K,  K1 ,  < 

(i<M2>-  IKJ  > 

66  t1?=M?M 
WO  112  =  1 

DO  30  III  1 ,  IM,  :< 

DO  30  F  1 ,  v 
K1=K-s  6 
til --Ml 
IKM2)  0 

DO  20  J  K,  Kl.  < 

IKM2>-  IKIIx)  H  <.j>  iiKiti) 
xii  h*=M*it 
<0  MX=MX 1 1 

II  <N  1)  9!5,  OS, 

•>:i  tli=3 

no  in  on 

96  til --9 

90  DO  91  J  1,111 

91  C<J)-|KJ) 

Kl  1UKN 
HID 
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SUBROUTINE  kt  III' 

c 

c-  -  ■  ?  JANUARY  <RFVISFD> 

C 

C  -  •  THIS  SUFTROU I  /  NF.  kl'flDS  IN  lilt:  I'llkllllF  Tl:l<s  I  uR  1  MCI  I  OKI  I  CHI.  SURFACE. 

C 

COMMON  kl'lV < 4 ,  ?!itJ > ,  DhVO,,  IRKX50),  DV < 50 ) ,  NBI  K,  Nk<24> 

DIMENSION  KH <  4  >,  VR <41 
C 

2  FORMAT  <4  <14,  1:1  4.  <0> 

4  FORMA  1  (1 110,  4X2  <IIOl‘  I  ICllI  sll|<Flii:h  NUMol  K  ,  l.t//<b<  IS,  iH>,  E14.  7))> 

6  FORMAT  <  tH  ) 

C 

C  •  ■  I  III  TIAI  1/F  |||F  OV  likkAV. 

C: 

IK  =  O 

:m  ik  =  ik  u 

DO  ??.  I  «  1,  26 
12  DVU>  «  0.0 

DO  33  I  30,  tio 
3<  DV<  I  >  0.0 

C 

C  -  kino  Dfnn  Fuk  lilkliUJlMK  (irillfll  sURHu.F  H<h||  ClikD  kFlfiDF.R. 

c 

14  Rt  AD  ?,  <  Kl'{<  1  >,  VB<  I  ),  1  1,4  ) 

DO  38  l  •  1, 4 
KK  =  KHCO 
IF  <  KK  >  3o ,  <k.  <8 
"iti  DV < Kis )  VB<  I  ) 

II  <  KK  F,n  )  <8,4  0,4  0 
’<8  CUNTINUF 
00  TO  34 
C 

C  •  81  CiRI-  01  I  I  till  MU'lin  I-  I'llRlIilt-  1 1  KS  nf|  dj'.i- 

C 

40  Clll  L  FMK1IE  <  4IIMok,  JK.  DV.  100  > 

42  II  <  IIF,  i  fill  <4l)M0k>  .Ns-  1  )  ofi  i  Cl  4r' 

II  <  DV (.sO  .IT  DV<29)  >  mm  ID  <u 

III  TURN 
C 

t  -  kkINf  FiikilMl- ll-ks  Ul  ok  i  1 1 ;  1 1 1  silkl-Tn 
C 

FAIRY  I  isf 
NFDV  2  DV<2(0 
MOV  "  DV<2>0 

f'kIMI  4,  NFdV,  •  l,DVk1>.I  I,  MOV  > 

FfcINI  G 
RF  fUkH 
1  ND 
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THIS  SUBROUTINE  MILL  ACCEPT  INPUT  FROM  1  HE  TYPEWRITER  FOR 
RETRACING  A  BEAM  THROUGH  THE  PRESENT  OPTICAL  SYSTEM.  IN 
ADDITION  THE  INPUT  BERM  MILL  BE  RELOADED. 

COMMON  RAY(4, 250), DEV(5,  250),  IRF(250), DV(50), NBLK, NR(24> 

2  FORMAT  (  1H@, 32HCHAN0ES  FOR  OPTICAL  SURFACE  NO.  ,127 
115HXX  :  XXXXXXXX  :  XXX,  12X,  1H<?  ) 

4  FORMAT  <  1H@,  14HINP1JT  FORMAT  N7 
133HXX  .  (VARIABLE  LOCAriON)  7 

230H : XXXXXXXX : XX  (VARIABLE  VALUE)  7 
325HX  (1-MIS TAKE)  ) 

8  FORMAT  (  12,  F.12.  7,  II  ) 

10  FORMAT  (  1H0, 26H INPUT  MISTAKE  RETYPE  INPUT  7 
115HXX : XXXXXXXX : XXX,  12X,  1H0  ) 


RELOAD  INPUT  BEAM. 

CALL  BEAM  IN 
IK  =  0 
RUN  =  8.  0 

IF  (  SSHTCHF  <6  >  .  EQ.  2  )  HR  HE  (53,4  ) 
INPUT  CHANGES  THRU  TYPEHRITER  PER  SURFACE. 


18  IK  =  IK  +  1 

CALL  F  READ  (  4HS I  OR,  IK, DV, 1 00  ) 

28  IF'  (  IRS  rATF(4MST0R)  .  NE.  1  )  GO  (0  20 

IF  (  IK  .  GE.  1  .  AND.  RUN  .  GT.  0.  0  >  DV(,'7) 

NS  =  DV(26) 

MRITE  (59,2  )  NS 
22  READ  (58,3  )  J,  A,  MISTAKE 

IF  (  MIS  I AKE  .  NE  1  )  GO  TO  24 
WRITE  (59,10) 

GO  TO  22 

24  IF  (  J  .  EQ.  50  )  00  TO  26 

IF  (  J  .  EQ.  27  )  RUN  =  A 

DV( J)  =  A 
00  TO  22 

RESTORE  UPDA1LD  SURFACE  PtlRAMt-  I LRS 

26  CALL  FWR1TE  (  4llb  TOR,  I !(,  DV,  1 00  ) 

20  IE  <  I WS  I  ATF  ( 4HS I  OR )  .  Nt-  1  )  GO  10  28 
IF  (  DV(76)  .IT  DV(29)  >  00  TU  18 
RETURN 
END 
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SUBROUTINE  SCAN  C  R  ) 

•  •  10  JANUARY  1S?4  CREVISl-D) 

- I'MIS  SUBROUTINE  DEVELOPS  (HE  DATA  NECESSARY  10  PLOT  ft  KNIFE  EDGE 

---  SCAN  OF  ft  B Eft II  ftr  THE  OUTPUT  REFERENCE  PLANE  OF  A 
■  '  PARTICULAR  OPTICAL  SURFACE. 

COMMON  RfiV < 4 >  2S0),  DEVC5/  2S0)/  1RFC250),  DVC50),  N6LK/  NRC24) 
COMMON/'Pl.TS/2IBUFCC.000)..  AC8)/  PC3),  XCSV,  YCS)/  I  NEC2)/  LINX/  LABLC5/  8) 
1  ,  XPC20O0)/  XMftX,  XMIi'L  DX/  YMAX/  YM1N/  DY 

DIMENSION  CC230,  2) 

EQUIVALENCE  CDFV,  C> 

2  FORMAT  <?0HKNIE  H  EDGE  SCAN  ) 

4  FORMAT  C20HPFR  CENT  ENERGY  ) 

6  FORMAT  C  1H@, 30HTYPE  1  FOR  YES  ,  lYPH  0  FOR  MO  / 

1  42HRKPI  OT  KNIFE:  EDGE  SCAN  Willi  NEW  DIMENSION!  ) 

8  FORMAT  <•  II  ) 

10  FORMAT  <  1H0/26HIF  MISTAKE  OH  INPUT  TYPE  1  X 

1  42HNEW  DIMINSloN  S  :  XXXXXXXX :  XX  ,  MiSltiKE  \  X  ,  IX/ 1HC*  ) 

12  FORMAT  <  E12  f,  It  1 

14  FORMAT  <  1H0»,  2<>IUNrUT  MIS  I AKE  kEIVPE  INPUT  > 

••  INITIALIZE  hie.  ARAl’II  PTiRAME  i  EkS  Mir  THE  KNIFE  edge  SCAN  PLOT. 


ENCODE  <  26/ 2.  I  HI51  <t,  1  >  ) 
ENCODE  <  20/ 4/ 1  AMI  <1, 8 1  1 
PCD  -  1.  0.1  *  PC2) 

YC?)  0.  SO  t  YC8) 

I  NEC1  1  0  LI  NFC  21 


PCI) 

YC?) 

I NE  Cl) 
YMIN  -■ 
YMAX  = 
DY 


0. 

1.00. 

20. 

0VC2O 
-  OVC?  O 
DVC>  :)/•!  0 


■  •  CALCINATE  Till  T  O I  I'll  NUIIIIEK  OF  kftVs. 

NRY  =  0 

DO  22  J  ■'  1/  NP.l  K 
22  NRY  -  NRY  +  NRCJ) 

YTOTHL  -  NRY 

-  FROM  THE  SCAN  INPKIMFNl  Clll  CUI  M  IF  IMF.  NUIIBElK  OF  POINTS  FOR  THE 
-  -  Pt  OT  AND  INI  I  I  III  I2F.  THE  X  AND  Y  IuSI1iiTN'> 


2  {  DSCAN  -  2  0  *  XMAX/2‘i0.  0 
N  -  2‘o0 

DO  24  J  -  t, N 
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ANULT  ~  J  -  1 
CCJ, 1>  -  XMIN  1  AI1UI.T=«D3CmN 
24  CCJ,  2)  =  0.  G 
C 

C-  SCAN  THROUGH  KOCH  BLOCK  OK  RAYS  AND  luiftl.  THK  NUMBER  OK  RAYS 
C—  FOR  EACH  X  POSITION. 

C 

DO  23  J  1,  NBI  K 
NRY  =  NR< J) 

IK  <  NRY  .EG.  0  >  00  TO  28 

CALL  FRK.flD  <  4HRfiYF>  J,  RAY,  8*NRY  > 

26  IK  C  IRSIRTI-"C4HRAYE  >  .HE.  1  >  GO  TO  26 
DO  28  I  -  1, N 
DO  28  K  «  1, NRY 

IK  <  RAYC1,  K>  .it-.  CCI,1>  )  CO,2>  =  CO, 2)  t  1.0 
28  CONTINUE 
C 

C  *  -  CALCULATE  THE  PER  CENT  OF  ENERGY  AS  THE  I  IGHl  BEAM  Ii  SCANNED 

C 

DO  30  J  =•=  1,  N 

30  CCJ,  2)  100.  *t:CJ,  ?>/YlOllTI 

CALL  F WRITE  <  4HXARK,  2*N  > 

32  IF  C  IW3 I  ATE  <4IIX'IRt  )  .  NK  1  T  UO  iO  32 
CALL  FWRlfE  <  4UYGRF, 1, CC1, 2), 2*N  > 

34  IF  <  I WS I MTK C  41 IVURK  >  .  ME.  .1  >  00  TO  34 
IBLK  =  NBLK 
INRY  -  NR  C 1 ) 

NBLK  =  1 
NR  C 1 )  =  N 
CALL  GRAPHS  <  R  > 

NBLK  =  JBLK 
NRCl)  -  INRY 
C 

C  -  ■  OPTION  TO  RM-I.ill-  GRAPH  MI  III  NEW  si:Al.E 

C 

WRITE  CSV,  6  > 

READ  (U8, 8)  I 

IF  C  I  .  EQ.  0  >  00  10  4o 

WRITE  <‘.<4,  10  >  • 

36  READ  CS3,  12)  D I II,  III  ST  AKE 

IF  <  .MISlAKK  ED.  A  )  00  lO  33 
WRITE  CSV, 14) 

GO  TO  36 

33  XMAX  -•  DIM 

XM I N  -  DIM 

DX  “  DIM/4.  A 

GO  TO  ?.< 

48  RETURN 
END 


I 
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c 

c— 

c - 

c 


THIS  SUBROUTINE  CRl.CULt'llES  THE  SINF.  AND  COSINE  FUNCTION  OK 
RNOLE  fl  <  IN  RADIANS  ).  THE  OU'iPUT  IS  IN  THE  B  ARRAY. 

DIMENSION  B(«) 

T  =  0 

B<1)  =  SINO  ) 

B<2>  =  COS  C  T > 

RETURN 

END 
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SUBROUTINE  TRACE  <  NkY, flZERO  ) 
•-  24  NOVEMBER  i9?.i  (REVISED) 


-  THIS  SUBROU f I NE  TRUCES  fl  ROY  FROM  RN  INPUT  RFFEktNCE  PLANE  THROUGH 

-  fl  REFLECTING  OR  REFRfiCtING  PLANE  SURf  ITCH  TO  RN  OUTPUT  REFERENCE 

-  PLANE. 


-  THE  C  RRKRY  IS  USED  FOR  INIERMFDIATE  CRI  CULfll IONS. 

-  RIL  RNGI  E  DATA  MUSI  BE  IN  RADIANS. 


COMMON  RAYC4,  250),  DEV (3,  250),  IRFOSO),  DVC30),  NBI.K,  NRC24) 
DIMENSION  C ( 40 ) 

DO  60  J  •••  1,  NRY 
IF  <  1RFCJ)  .  KQ.  1.  )  Gu  lu  60 
IF  <  DVC22)  .  FO  0.  0  )  GO  TO  35 
DO  30  I  =  1, 3 
30  DVCI+4)  -  DI-VC  I,  J) 

35  CONTINUE 

DO  40  I  -  1, 40 
40  CCI  )  =  0.  0 


CRLCUlflTE  DIkT.CIION  COSINES  OF  INPUT  kflY. 


T 1  =  TONI*  C  RAYC3, J)  )  $  T2  -  TAME  C  kAY(4,  J)  ) 

CC30)  =  COS  C  mON  C  SORT  <  TT+T1  f  T/f’I  V  >  >  ) 

CC28)  =  V1-FCO0)  $  C < 29 >  *  T2rCO0) 

CC31)  =  RRY  (1)  J)  $  C<32)  =  RAYO.J)  *  C<<2)  =  -DVC2) 

C<5)  “  DV<1) 

DO  42  I  “•  1,  2 

42  CALL  SI  NCOS  j.  C  DVCI  F4  ) ,  CC2  F I -1  )  ) 


1=1,  CflLCAl  A I K  POINT  Of  kEf  I  ELI  ION  uk  REncrtCi  ION. 

1=2,  CflLCAl  Alh  POJNl  OE  RflY  ON  OU I  PUT  PLANE  AND  ITS  DIktCTlON. 


DO  50  i=-:,2 

CALL  SINC0S1  C  CCS)fDV(3), C<6)  > 


IF  C 

DV<3>  .  f  Q. 

O 

0  )  CC?) 

=  A/FRO+C  <  T ) 

CC10) 

=  C<4) 

* 

CC13)  = 

C(3)fC(?) 

T 

C<)6>  = 

CC3)fC(6> 

C(ll) 

Q 

© 

n 

* 

C  C  I  4  )  -- 

-CC6) 

t 

CCI 2)  = 

CC7) 

CC12) 

=  C(3> 

* 

CO  5)  = 

-CC4)t<CCF) 

* 

COW)  -- 

-*CC4)>fC(6) 

IF  < 

I .  fc  C  2  )  OA 

TO  52 

CC19) 

=  C  <  4  ) 

t 

Cl??)  ~- 

CC  OfCCl) 

* 

C(?3>  ~ 

-C(3)fC(2) 

C<20> 

=  0  0 

1 

C  C  2  3  )  = 

CC2) 

t 

CC  V6>  - 

Ci.1) 

•;c2i> 

=  CC3) 

'T 

CC24)  = 

CC*)fC(1  > 

1 

CO'.')  - 

CC  4  ) i'C C 2  ) 

00  TO 

54 

CALL 

SINC031  < 

CCS), CCA)  ) 

CO) 

=  A/KRO^CO) 

CC19) 

=  J  0 

i- 

0(2?)  ~ 

0  0 

t 

COS)  - 

0.  0 
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C<28>  =  8.  0  $  C<23>  =  C<8>  *  C<26> 

C<2±3  =0.0  *  C<24)  =  -C<9>  *  C<??> 

54  CALL.  MATVEC  <  C<3 8),  C<28>,  C<28>,  I-i  > 

CALL  MATVEC  <  CCIO),  C<31>,  C<31  >,  I--i  > 

CALL  MATVEC  <  C<19>, DV<?>  ,  C<34>,0  ) 

DO  56  L  -  31>  33 
36  C<L>  =  C<L)  -C<L.+3> 

IF  <  I  .  LQ.  2  >  CCCO  =  C<3i)-0VCO 
T1  =  C<2S)/C<30> 

T2  =  C<29>,  J<30> 

C<3±>  =  C<31)  -  T1*C<33> 

C<32)  =  CG32>  --  T2*C<3:0 
C<33>  =  0.  8 
IF  <  I.  EG!.  2  )  GO  TO  50 
CC28)  =  DV<31>t=C<28> 

C<23>  =  0V<11.)^C<29> 

CG30J  =  •  DVCI.0>*SQI<IF  <  1.  0  DV<11  >*DV<±i>*<  3. 
CG5)  =  DVG'0>):nCriSF  (  DV<±i>*COS  <  C<5>  >  > 
50  CONTINUE 


C 

C- 

C 


•-  OUTPUT  RAY  RELATIONSHIP. 


RAYCLj J) 
RAY<2« J) 
RAY<3, J) 
RAY<4# J) 
60  CONTINUE 
RETURN 
END 


C<3 3  ) 

C<32> 

RIAN  <  11  ) 
A  I  AN  <  12  > 
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SUBROUTINE  TRANSFER 
C 

C  24  NOVEMBER  1972  <RKVISED> 

C 

C  THIS  SUBROUTINE  MILL  TRANSFER  THE  RHV  DAI  A  GENERATED  BY  THE 
C  RAYTRACE  PROGRAM  FROM  1  HE  RAYF  FILE  TO  1  HE  BEAM  FILE.  THIS  WILL 

C  ONLY  BE  NECESSARY  IF  IT  IS  DESIRED  TO  SAVE  THE  DATA  FOR  LATER  USE. 

C 

COMMON  RflY<4>  250>,  DEV<5,  250>,  IRF<250>,  DV<50>,  NBl.K,  NR<24> 

CALL  F WRITE  <  4IIMRYF,  3  ,  NBI.K,  25  > 

52  IF  <  IWSTATF <4HNRYF)  .  NE.  1  )  GO  TO  52 
DO  60  K  =  1,  NBl-K 
NRY  =  NR<K) 

IF  <  NRY  .  EQ.  0  >  GO  TO  60 
CALL  FRHflD  <  4HRRYF/  K>  RAY,  8*NRY  > 

34  IF  <  IRS  I  RTF<4HRAVF  >  .  NE.  3.  >  GO  10  54 
CALL  FWRI'fE  <  4HBEAM,  K,  RAY/  8+NRV  ) 

56  IF  C  I  US  I'HTF  <4  HUE  AMD  .  NE.  1  )  GO  10  56 
60  CONTINUE 
RETURN 
END 
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SUBROUTINE  WRITE 

-  9  NOVEMBER  19?.? 

-  THIS  SUBROUTINE  Will.  WRITE  THE  OUIPU1  BEAM  PARAMETERS. 

COMMON  RAY<4,  250>,  BF.V<5,  250>,  IkF<250>,  DV<50>,  NBl.K,  NR<24> 

4  FORMH7  <  1H",  4X11H0UTPUT  BERM,  1H:  > 

6  FORMAT  <15,  111),  4F.20.  ?> 

8  FORMAT  <  1H0,  4X1  5HBI.0CK  NUMBER  =  ,  12, bXJ  <HBI  uCK  SI2E  =  ,12//  > 
PRINT  4 

DO  26  I  =  1,  NBIK 
NRV  -  MR< I ) 

IF  <  NRV  .  EO.  0  >  1,0  TO  26 
CALL  FRFAD  <  4HRAYI-,  I,  RAV,  8+NRY  > 

22  IF  <  I  RSI  ATF<4HR0VI- )  .  NF.  1  )  GO  10  22 
PRINT  8,  I, NRV 
DO  24  J  =•-  1,  NRV 

24  PRINT  6,  <  J, <  RAY<K, J ) ,  K“i, 4  )  > 

26  CONTINUE 
RETURN 
END 
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